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Abstract 

We continue our study of the retarded Green's function of the universal 
fermionic supersymmetry current ( "supercurrent" ) for the most general 
class of d = 3 N = 2 SCFTs with D = 10 or D = 11 supergravity du- 
als by studying the propagation of the Dirac gravitino in the electrically 
charged AdS-Reissner-Nordstrom black-brane background of N = 2 min- 
imal gauged supergravity in D = 4. We expand upon results presented in 
a companion paper, including the absence of a Fermi surface and the ap- 
pearance of a soft power-law gap at zero temperature. We also present the 
analytic solution of the gravitino equation in the AdS2 x M 2 background 
which arises as the near-horizon limit at zero temperature. In addition 
we determine the quasinormal mode spectrum. 
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1 Introduction 

All N = 2 SCFTs in d = 3 dimensions have a universal current supermultiplet 
consisting of the energy- momentum tensor, the abelian global i?-symmetry current, 
and the fermionic supersymmetry current or "supercurrent" . We will be interested in 
the general class of such SCFTs that have weakly coupled supergravity duals in either 
D = 10 or D = 11 dimensions. We will carry out a comprehensive investigation of 
the Green's function of the supercurrent for this universal class when held at finite 
temperature and non-zero chemical potential with respect to the global .R-symmetry. 

A significant motivation for this paper, which technically clarifies and conceptually 
extends the companion paper [I] , is the extensive work on fermion spectral functions 
that has been carried out using gauge-gravity duality starting with [2jj5] and then 
further developed in [6 56 . Remarkably, all of these works have, essentially, been in 



a bottom-up context, with the exception of [23] which studied the spectral functions 
in the context of probe-brane p-wave superconductors. It is therefore of considerable 
interest to carry out more detailed top-down analyses to see which features can ac- 
tually be realised in a string/M-theory setting. An appealing feature of the analysis 
that we carry out is that it is applicable to a universal class of SCFTs. Note that 
after the appearance of [l] and as this paper was being written up, a related study 



appeared in 57 



Consider, then, the most general AdS& x M 6 and AdS± x M 7 solutions of D = 10 
and D — 11 supergravity, respectively, that are dual to N = 2 SCFTs in d — 3. 
There is very strong evidence that there is always a consistent KK truncation of the 
D = 10 or D = 11 supergravity theory on M 6 or M 7 , respectively, to minimal N = 2 



gauged supergravity in D = 4 58 . This means, by definition, that any solution of the 
latter theory can be uplifted on any of the Mq or M7 to obtain an exact solution of 
D = 10 or D = 11 supergravity, respectively. Indeed for the infinite class of D = 11 
solutions when M7 is a seven-dimensional Sasaki-Einstein manifold, SE7, and also 
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for another general class of D = 11 solutions with magnetic four-form flux, this was 
proven at the level of the bosonic fields in [58] and, furthermore, consistency with the 
supersymmetry variations was shown. For the SE7 case it has also been shown that 



the truncation is consistent in the fermionic fields, at the quadratic level, in 59 , 60 



The bosonic fields of minimal N = 2 D = 4 gauged supergravity consist of the 
metric and an abelian gauge field, which are dual to the energy-momentum tensor 
and the R-symmetry current in the N = 2 d = 3 SCFT, respectively. The D = 4 
supergravity also has two Majorana gravitini, or equivalently a single Dirac gravitino, 
which is dual to the supercurrent in the SCFT. An important feature of minimal 
N = 2 D = A gauged supergravity is that the bosonic Lagrangian is simply Einstein- 
Maxwell theory with a negative cosmo logical constant. In particular, the electrically 
charged AdS-Reissner- Nordstrom (AdS-RN) black-brane solution solves the equations 
of motion. Thus, the results of [58] imply that this solution is relevant as a dual 
description for a/Qof the d = 3 N = 2 SCFTs with gravity duals at finite temperature 
T and non-zero chemical potential // with respect to the R-symmetry. The strategy 
for calculating the Green's function for the supercurrent is to solve the linearised 
gravitino equation in the AdS-RN background and then examine the asymptotic 
behaviour at the boundary. To calculate the retarded Green's function, which will be 
our main objective, we need to impose ingoing boundary conditions on the gravitino 
solutions at the black-brane horizon. 

The Green's function that we obtain is valid for the entire class of N = 2 SCFTs 
at least for large values of the dimensionless temperature T/fi. However, as T/fi is 
lowered it is possible that the field theory undergoes one or more phase transitions. 
The corresponding phases will be described by new black-brane solutions which will 
involve KK modes lying outside the truncation to minimal gauged supergravity and 
hence will depend on the details of the specific SCFT that is being considered. For 
example, in the SEj class of SCFTs it is known that there are superfluid instabilities 
which spontaneously break the R-symmetry 61 . An analysis for the case of S 7 can 
be found in 62 and for another class of SCFTs in 63 . It would certainly be very 
interesting to know the precise nature of all of the phase transitions all the way down 
to their ultimate zero temperature ground states, even for specific examples of M 6 
or M7 (for some discussion see (62])- One could then aim to extend the calculations 
of this paper and [I] by solving the gravitino equation in the thermodynamically 
preferred black-brane background for any given temperature. On the other hand, it 



1 Of course the electrically charged AdS-RN black-brane can also be relevant for certain non- 



supcrsymmetric CFTs. A prominent example is the "skew-whiffed" AdS^ x SE7 solutions 61 
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is possible that for at least some N = 2 SCFTs the AdS-RN black-brane is valid all 
the way down to zero temperatures. 

The paper is structured as follows. In section [2] we begin with some general 
comments concerning the Green's function of the supercurrent including a discussion 
of discrete symmetries. In section [3] we introduce N = 2 D = A gauged supergravity 
while in section [4] we derive a convenient form of the gravitino equations of motion. 
In section [5] we analyse the asymptotic behaviour of the gravitino equations at the 
AdS boundary and derive a useful formula for the Green's function. 

We will calculate the retarded Green's function and hence we must impose ingo- 
ing boundary conditions on the gravitino equations at the black-brane horizon. The 
details of how to do this for non-zero and zero temperature are different and so are 
treated separately in sections [6] and [7j respectively. In both cases we solve the equa- 
tions numerically. We will see that the spectral function, as a function of frequency 
u and wave-number k = |k| has a "phonino" pole at u> + fi = and k = 0. This is a 
collective fermionic excitation that any supersymmetric medium at finite temperature 
has. This long-wavelength mode is important for the hydrodynamical description and 
has been studied, for zero chemical potential, in 64 66 . Furthermore the dispersion 
relation of the pole was obtained holographically in 67 . Our work extends these 
investigations to non-zero chemical potential. We will see that the pole persists at 
T = and we will also extract the dispersion relation. It is worth emphasising that 
in our conventions u = corresponds to the Fermi energy and hence this phonino 
pole is not gapless^] when /z ^ 0. 

Another important feature of the Green's function for T ^ is a depletion of 
spectral weight near u = 0. At T = the spectral weight vanishes at u = with, 
for each k, a power-law behaviour. This latter feature is controlled by the locally 
quantum critical theory dual to the AdS2 X K 2 background which arises as the near- 
horizon limit of the AdS-RN solution at T = 0. To see this in detail, in section [8] we 
first obtain the exact solution for the gravitino equations in the AdS2 xl 2 background. 
Then in section [9] we show how the non-analytic structure of the Green's function 
for oo ~ 0, as well as the power-law scaling, is governed by the exact results in the 
AdS 2 x M 2 background. 



In section [10] we calculate the spectrum of quasinormal modes for T ^ finding 
results consistent with some of the main features of the spectral functions. 



We conclude in section 11, summarising some of our main results. It is worth 
emphasising that we find no Fermi surface for the supercurrent correlation function 



2 This is reminiscent of the ferromagnetic magnons discussed in 
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in our supergravity calculation. We suspect that in the top-down models that we are 
considering a Fermi surface may be seen in other fermion correlators and / or at higher 
loop order. Another interesting feature is that we will see the conformal dimension 
of the operators in the one-dimensional conformal field theory in the far IR, dual to 
the AdS2 region, are all real. This means that the log-oscillatory behaviour see in 
the bottom-up models of [3|[5] is absent. It is also a manifestation of the fixed charge 
over mass ratio dictated by the string/M-theory compactifications we consider. In 
physical terms this charge to mass ratio is never high enough to make pair production 
of charged Fermions energetically favourable IB"]. This leads us to speculate that the 



system displays no instability towards forming an electron star 18 35 53 , a state 
which can be viewed as the result of occupying the corresponding modes in the bulk. 

Some technical material is presented in four appendices. Appendix [A] contains 
our spinor conventions. In appendix [B] we discuss the positivity properties of the 
spectral function. Appendix [C] contains a calculation for the asymptotic behaviour 
of more general Rarita-Schwinger equations in arbitrary spacetime dimensions, but 
propagating in pure AdS space. Although not strictly needed for the main text we 



found the calculations enlightening. In appendix D we summarise some aspects of the 
numerical methods that we employed to obtain our results. Appendix|E]contains some 
technical material that we use in calculating the Green's function in the AdS2 x M. 2 
background. 



2 The supercurrent correlation function 
2.1 The form of the correlation function 

Let S a be the conserved supercurrent operator of the dual d = 3 N = 2 SCFT in IR 1,2 . 
It is a complex vector-spinor: the index a, with a = t,x, y, is the vectorial index and 
we have suppressed the spinor index. Since the supersymmetric theory is conformal, 
the supercurrent is gamma-traceless. We aim to calculate the retarded correlation 
function G a p(p) = ^S a (p)Sfi(0)) R at temperature T and chemical potential \i with 
respect to the abelian i?-symmetry. We will do this by exploiting the fact that the 
expectation value of the supercurrent in the presence of a vector-spinor source a a , at 
linearised order in the source, is given by 

(S a )=iG aP a p . (2.1) 
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The gamma-tracelessness and conservation of the supercurrent implies that 



7 Q (5 a ) = 0, p a (S a )=0, (2.2) 

where 7° are d = 3 gamma-matrices. Our conventions for d = 3 and D = 4 spinors 
can be found in appendix A. Since we are considering the SCFT at finite /x, which 
can be viewed as weakly gauging the i?-symmetry, we have 

p a = {u,k), u = u + fi. (2.3) 

Note that u = corresponds to the Fermi energy. Since the source couples to the 
supercurrent by a term in the Lagrangian of the form J d 3 x(S a a a + a a S a ), the source 



in (2.1) can be taken to satisfy 

1 a a a = 0, Sa a = (Si - §7«7 /3 ) Ppe , (2.4) 



where the second equation in (2.4) arises from the weak gauging of the supersym- 
metry. Of course the supercurrent itself, and hence its expectation value, is gauge 
invariant. In general, since S a has two independent components, G a p will have four. 

It will be convenient to extract the four components of G a p by introducing a basis 
of 3d vector-spinors ea\ i — 1, 2, satisfying 7 a ea = p a e£) = and the orthogonality 
condition e« e^ a = — 2p 2 e u . We can then write 

G a ? = tyegef , (2.5) 

where the four functions of p a are the four independent components of G a p. Note 
that the mass dimensions of G a p{p) is 2 and hence tij is dimensionless. We also ob- 
serve that the expectation value, (S a (p)), and the source, a a (p), have mass dimensions 
— 1/2 and —5/2, respectively. 

An N = 2 SCFT in R 1 ' 2 at finite temperature and chemical potential will preserve 
0(2) rotations (and reflections) in the x — y plane and this strongly constrains the 

as we now show. We first rotate so that 

p a = (u,k,0), (2.6) 

with k = |k| and take the following explicit basis 

e£\u>,k) = (k, -w,0) a n-(0,0, k-u) a m, 
e£\u, k) = - (k, -u>, 0) a m + (0, 0, Q + k) a n , 
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where m, n are two two-component spinors satisfying / y tx m = m and ^ tx n = —n. We 
now consider the d = 3 parity operator that takes (t,x,y) — > (t,x,—y). Using the 
conventions described in appendix A we have 

C W(w, k) ^(ii, fc) = -e«(w, fc) , 

e^(u>, fc) Vj-feffa k) = e «(a>, fc) , (2.8) 

where "Pq,' 3 = diag(l, 1, —1). This induces the action £i 2 — >■ — £12 and £ 2 i — >• —£21 at 
the same value of (a), fc) and hence we can conclude that 

£ 12 = t 21 = 0. (2.9) 

We next consider the 50(2) rotation that takes (t,x,y) — > (t,—x,—y). A short 
calculation shows that 

e%\u;,k)^e a l \u;,k), (2.10) 

and hence 

t 72 (u>,k) = t 11 (u,-k). (2.11) 

We find that all of the information of the retarded Green's function is contained in 
tn (a), k), if it is formally viewed as a function of both positive and negative k. 

Thus our objective is to calculate tn(u),k) and more specifically the spectral 
function A(Cj, k) which we define to be the imaginary part: 

A(u, k) = lmt u (uj, k) . (2.12) 

In appendix [B| we show that A(u, k) > 0. 

Before concluding this subsection, it is interesting to relate the above discussion to 
an alternative way of characterising the four independent components of the general 
Green's function. We first introduce the projector onto transverse gamma-traceless 
parts of vector spinors 

PJ =&c?-\ {la - P ~f)l"-^ (3p« - 7«|0 / • (2-13) 
Then the constraints ( |2.2 ) imply that we can write G a p(p) in the general form 



G aP {p) = P Q 7 [t l5pl p + 2{t lP 2 + t 2 p)Vjs] P 
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P 



£ a /3 7 7 7 + h [eap^f + {p 2 r] a p - p a pp)] (2.14) 
+ h [(p 2 i] al3 - p a pp)f + p 2 e a / 3l p' y ] , 
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where t a/3 ^ is totally symmetric, and satisfies p a t a /3^ = and t a a p = 0, and hence 
has two independent components. In general the four independent quantities t a p 7 , 
t\ and £2 are arbitrary functions of p a . In the case of interest, where we have 0(2) 
invariance, we can use the above results to conclude that 

tafr = h = 0, t 2 = , 1 u tn(u),k) . (2.15) 



This easily follows by sandwiching (2.14 ) between fh, n and m, n. If there is full d = 3 
Lorentz invariance we would have t Q( g 7 = t\ = and £2 = hip 2 )- Finally, we note 
that we also have ei'e^E^ = — 4p 2 P a /5 . 

2.2 Phonino 

Placing the N = 2 SCFT at finite temperature breaks supersymmetry and leads 
to a fermionic Goldstino mode. When [i = it was pointed out that this mode 
appears as a pole in the spectral function for the supercurrent correlator located at 



uj = 0, k = 64 66 . This mode is important in the hydrodynamical description 
of the system 64 , 65 and in this context is known as the "phonino" . Recall that to 
construct the hydrodynamics one needs to identify a set of variables whose thermal 
expectation values distinguish equilibrium states of the theory. These consist of the 
charge densities for all conserved currents. In the present setup there are the usual 
energy and momentum densities associated with the energy-momentum tensor, there 
is the R-charge density associated with the R-symmetry current, and there is also the 
supercharge density associated with the supercurrent. The hydrodynamical equations 
are obtained by postulating constitutive relations for the spatial parts of the currents 
in terms of the charge densities and then imposing current conservation. By analysing 



the constitutive relations when /j, = it was shown in 65 that in addition to normal 
sound, there is also, for small k, a weakly damped propagating fermionic collective 
excitation. This phonino or supersound is somewhat analogous to second sound in a 
superfluid. The phonino was first studied from a holographic point of view at finite 



T and [i = in 67 for SCFTs in d = 4. In particular, the diffusion constant for the 
dispersion relation for the phonino pole was calculated. 

Here we will be considering finite temperature and fi 7^ 0. In view of our comments 
above, we expect that switching on // will move the phonino pole to 00 + ^1 = 0, k = 0, 
and our results will confirm this expectation. Recall that 00 = corresponds to the 
Fermi energy and hence the phonino is certainly not a gapless excitation. From our 
numerical results we will also be able to extract the dispersion relation for the phonino 
pole. 



S 



3 N = 2 minimal gauged supergravity in D = 4 



The field content of minimal N = 2 gauged supergravity in D = 4 69 70 consists of a 



metric, a gauge field A and a single Dirac gravitino ip/j, (or equivalently two Majorana 
gravitini). The action is given by 



S 



d A x\/—g 
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(3.1) 



where T = dA and we have not explicitly written the four-Fermi terms, which will 
not be needed. Here 



D = V - iqA, 



-m 



(3.2) 



where V is the Levi-Civita connection. Observe the presence of Pauli terms which 
couple the gravitino non-minimally to the field strength of the gauge fielcj^j It is also 
useful to record the supersymmetry transformations about a bosonic configuration. 
The only non-trivial transformation is for the gravitino and is given by 



777 ? 

8%= (^--r M + -^ p r^r M ] 



(3.3) 



where e is an infinitesimal Dirac spinor. 

The AdS-RN black-brane solution solves the equations of motion and is given by 



ds 2 



-fdt 2 + 



dr 2 r 2 



f 



A = <pdt . 



with 



(3.4) 



9 / 2 \ 2 

<Y* * ly / 7™ \ <Y' 

' 14-11+ ( fJ 2..2\ , i)2 ,.2 ' J 



The horizon is located at r = r + and the temperature, T = /'(r + )/47r, is given by 

T = (3r + /£ 2 - £ 2 /j 2 /r + )/A7i . (3.6) 

Thus, when ft = \^3r + /£ 2 we have T = and 

(r — r + ) 2 (r 2 + 2rr + + 3rjJ 



/ 



(T = 0) . 



(3.7) 



3 In a bottom- up context the role of Pauli term couplings for spin-| fermions have been studied 



in some detail in 38 43 and 46 . 
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Furthermore, as r — > r + the T = black-brane solution approaches AdS 2 X K 2 
with the radius of the AdS2 given by L( 2 ) = ^/\/6. Indeed after introducing a new 
coordinate via r — r + = L 2 2 ^/z the solution approaches the exact solution 



ds 2 



.A 



V2z 



dt. 



(3.8) 



We will study the equation of motion of the gravitino (at the linearised level) in 
the AdS-RN background ([34]), ([375]) and also in the AdS 2 x 1R 2 background (Q. It 



will be convenient to use the local supersymmetry transformations (3.3) to work in 
the gauge 

D^ M = r^ M = 0, (3.9) 
in which the linearised equation for the gravitino is given by 



fj,l/ 



ij p + iXF^T^T^ = . 



(3.10) 



Here we have introduced A = 1 to highlight where the Pauli coupling enter in some 



places in the sequel. Note that (3.9) does not quite fix the gauge completely. There 



are residual gauge transformations of the form (3.3) where e satisfies 



{V^Dp - 2m)e = 



(3.11) 



4 Effective two-dimensional formulation of the grav- 
itino equation 



In solving the gravitino equation (3.10) in the AdS-RN background (3. 4), (3. 5) and the 
AdS 2 x IR 2 background (3.8), it will be convenient to, effectively, carry out a dimen- 



sional reduction on the two spatial directions x l to two space-time dimensions. This is 
advantageous since the Lorentz representations of Spin(l,l) are all one-dimensional. 
It is convenient to consider the general metric and gauge connection of the form 



ds 2 = ds\ + E 2 dx % dx % , 

ds\ = g mn dy m dy n , A = A rn dy m , 

where H = H(y), g mn = g mn (y) and A m = A m (y). We choose a frame 



(4.1) 



Hdx 1 



(4.2) 
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in which the spin connection has non-zero components w'„ = e l dm In H and u) m h = 
Cj m h- We then use the planar symmetry of the background to expand the frame 
components of the gravitino in plane waves: 

Mv m ^) = j 0^ H ^(y m ^\ (4.3) 

where, for later convenience, we have included the factor of H' 1 in the definition of 
Xfi- We now focus on the two-dimensional m components of xp, since the transverse 
Xi can be determined via the conditions T^ip^ = D^ip^ = 0. 



Substituting into the equation of motion for the gravitino (3.10) we find 



[b + iH-%-m- \i\F np T n ^ xra + (din In H)T n Xn + i\F^T n T AXp = , (4.4) 

where j£ = fcjP. We would now like to expand Xm m terms of irreducible Spin(l, 1) 
representations. To do this we first note that we can decompose the spinor com- 
ponents under the two-dimensional chirality operator = T tr and also = 
{l/k)T^)k = r' Vf . This is possible since these two operators commute (for more 
details on our spinor conventions in D = 4 see appendix [A]). Let us first decompose 
under writing 

Xm = X m + X m ' (4-5) 

with T^Xm = +Xrk an d ^ X%, = ~X m - It wm be important to note that x m an d 



are not mixed by (4.4). Next consider the two-dimensional lightcone frame 

ds 2 = e + e~, e ± = e"±e i . (4.6) 

We can then introduce a basis for Xm given by 

x v _ u ( 3 /2) e + <g> r] + + u (1/2) e + ® rf + ?i ( ~ 1/2) e" ® r] + + w ( ~ 3/2) e _ <g> rf , 

x p = v ^l 2 ) e + <g> p + + v {1/2) e + ® p~ + v { - 1/2) e~ <g> p + + v { - 3/2) e~ g> p~ , (4.7) 

where = r^p 3 *" = ip^ and hence the superscript onu,t) refers to the 

helicity of the Spin(l, 1) representation. The equations of motion for x v then reduce 
to 

p_ n (3/2) + Ku am + ((9+ ln i? ) w (-i/2) = o ; 

V + u^ + ^V 3 / 2 ) + (d + \nH)u {1 ^ = , 

P_ w (-V2) + ^' u (-3/2) + ((9 _ ln ff ) w (-l/2) = ) (4 ' 8) 

P +M (- 3 / 2 ) + i^" 1 / 2 ) + (cL mif)^ 1 / 2 ) = , 



11 



where we have introduced the covariant derivatives 



V±u {s) = (d± - iqA± + 2stt±)u 



0J±. 



(4.9) 



and 



K 
K' 



[m + ikH~ l - 2i\F ] 
[m - ikH- 1 - 2i\F ] 



(4.10) 



where J 7 = F Q e + A e _ . The equations of motion for \ p give rise to equations for v of 
the same form but with k — > —k, or equivalently K <H- i^'. 

We can make a similar reduction for the residual gauge transformation. Decom- 
posing the Fourier modes as 

d 2 k 



e(x m ,x) 



(27T) 



.^-1^(1/2)^+ + ^(-1/2)^- + z (l/2) p+ + 



-l e ik-x 



(4.11) 



with u> and z functions of (x m ,k), one finds that the constraint which the residual 
supersymmetry transformations must satisfy, given in (3.11), can be written 

V-wW + P«; ( - 1/2) = 0. 
V + w^ + Pw^ = 0. 

with 



(4.12) 



P = m + 



2H 



(4.13) 



The equations for z are of the same form but with k — > — k or, equivalently, P <-> P. 



The supersymmetry variations (3.3) imply that 

5 M (3/ 2 ) = (I? + -a + lnP ( ) W ( 1 / 2 ), 
^^/^ = (D + - 9+ InP) + (|m + i\F ) 

5u { - 1/2) = (X>_ - <9_ In if) w (1/2) + (§m - iAF ) w 

5 M M/2) = (p__ 9 _l ni/ ) w (-l/2). 



(-1/2) 



(4.14) 



with similar equations for 8v^ s \ 

In the following we will just focus on the r\ dependent equations, namely the 
equations for given in (4.8) and the residual supersymmetry transformations 
given in (4.14), since the corresponding p dependent equations involving v^ s ' can be 
obtained via k — )■ —k. 

We now write down the explicit gravitino equations in the AdS-RN black-brane 



background (3. 4), (3.5). We will assume a time dependence of the form e tuJt . Adopt- 
ing the obvious orthonormal frame e* = f 1 ^ 2 dt, e r = f~ l l 2 dr, e % = {rjt)dx % we find 
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that the gravitino equations (4.8) take the explicit form 

3/' 



fl/2 



m H h iXcp 

r 



u (l/2) + = 

r 



« w + 



./ 



/' 1 

— + - 

4/ r 



r_ i 

/ 4/ + r 



(-1/2) 



/V 2 
1 



_z(ur+#) 3f 
/ + 4/ 



fl/2 



fl/2 

iki 



m h i\<p 

r 

m H 

r 



m i\cf)' 

r 



wC-Va) + = 



(4.15) 

We reiterate that g = — m = i~ x and A = 1. Observe that these equations are 
invariant under — > (u^^)* and we will return to this point in section 



9.3 



The residual supersymmetry transformations acting on the given in (4.14) 
take the explicit form 



Sit 1 '** 



fl/2 


d r - 


i(uj + q(j 


i) /' 


2 


f 




fl/2 


d r - 


i(uj + 


>) , /' 
4/ 


2 


f 


fl/2 

2 


d r + 


i(uj + qq 
f 


i) , /' 
4/ 


fl/2 


d r + 


i(uj + qq 


i) /' 


2 


f 


4/ 



(1/2) 



u^ 1 / 2 ) + [m - zAd r 0] w^) , 
(1/2) + i [m + t\d r <P] , 



w 



w 



(-1/2) 



(4.16) 



where the residual supersymmetry parameters u^ 1 / 2 ) satisfy (4.12) which can be 
written as 



r+ f + 4/ 
_ ^(u; + g0) f 
/ + 4/ 



^(V2) + 



1 



fl/2 
1 



fl/2 



2m + 
2m — 



iki 
r 

iki 
r 



w 



(-1/2) 



w 



(1/2) 



0. 



0. (4.17) 



Note that in these expressions u only appears in the combination (a; + q<p) which 
approaches (u + i~ x [i) at the AdS± boundary. This shows that u> = is the Fermi 
energy. 



5 Asymptotics and the Green's function 



In the remainder of the paper we will solve the gravitino equations (4.15), taking into 
account the residual supersymmetry transformations (4.16), (4. 17). As we are inter- 
ested in obtaining the retarded Green's function, we will impose ingoing boundary 
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conditions at the black-brane horizon located at r = r + . The boundary conditions 
for T 7^ and T = are different and will be discussed separately in later sections. 
The behaviour of the solutions at the asymptotic AdS boundary, located at r — > oo, 
allow us to extract the source and the expectation value of the supercurrent in the 
dual SCFT. The imposition of the ingoing boundary conditions will mean that the 



expectation value is fixed by the source, as in (2.1), and this allows us to extract the 



retarded Green's function. In this section we explain this procedure in more detail. 
5.1 Asymptotic expansion 

In analysing^] the asymptotic behaviour of solutions as r — > oo, it will be very help- 



ful to return temporarily to the covariant gravitino equations given in (3.9), (3. 10). 
Furthermore, it will be helpful to consider the radial components, ipf, and the d = 3 
components, ipa, of the bulk gravitino field separately, where hats denote D = 4 
tangent space indices. It is also helpful to further separate each of these into posi- 
tive and negative chirality components with respect to T r , denoted by superscript ±, 
respectively, in the following. At the end of this section we will reconnect with the 
equations for the t and f components that we derived in the last section. 



We find that as r — > oo the asymptotic expansion of (3. 9), (3. 10) takes the 
schematic form 

if, = r" 1/2 ^_ 1/2 + r- 3/2 ^_ 3/2 + r- 5/2 ^_ 5/2 + r" 7/: V_ 7/2 + (r- 7/2 logr)0_ 7/2 + . . . 

(5.1) 

In particular we find that log terms only start to appear at the order given. In more 
detail, introducing ± superscripts corresponding to positive and negative chirality 
with respect to T r , respectively, at leading order we find 





-1/2,Q 


= o, 




-1/2.Q 


= 


i>\ 


f 

-1/2,* 


= 0, 


i>'. 


-l/2,r 


= o, 



(5.2) 



with (1 + T r )A& = and T a A& = 0. This data is enough to obtain the next order of 

We found it useful, as a warmup, to first study a more general Rarita-Schwinger equation in 
arbitrary dimensions propagating in pure AdS space. The analysis is presented in appendix in (|C|). 
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the expansion and we get 

^-3/2,a = "^H + \^&{P ■ A) , 

4-3/2,& = > 
f 

-3/2,r 

-3/2,f ~ ~~ 2 

At the next order, new data appears. We find that 



^-3/2,f = , 

r. 3/2 .r = -¥^(P-A) (5.3) 



^-5/2,& = 




^-5/2,fi = 


lf[ p& ( p .A)-p 2 A & ] 


^-5/2,r = 


-\tf{p-A), 


^-5/2,r = 


-ie(g-A), 



(5.4) 

where (1 - F f )B & = and 

<te = (4w+,0,0). (5.5) 
New data also appears at the next order. We find 

^-7/2,& = W v & + 2P&P)(P-A) - 2p 2 pA & ] , 

^-7/2,f = ^ ' 

VT 7/2 ,, = ^V(P-^), (5.6) 
where (1 — r r )_B f = and we have defined 

= c diag(2, 1,1), c = -\{r% + t^r + ) . (5.7) 
Note that Ca a = 0. Finally, for the log terms we get 

#W = ' 

^7/2,a = ¥ 3T ^(9 ■ A) - \C^T & T^ , 
#W = - 2 ^(5 • A) + CtoTffa , 

4T W = . (5.8) 

At this point we have shown that the expansion is determined by the negative 
chirality, gamma-traceless the positive chirality B& satisfying various constraints 
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arising from (3.9), which we will not write down explicitly, as well as the positive 
chirality spinor B? . To proceed it is now very helpful to carry out a similar expansion 
of the residual gauge transformations. In particular, this will reveal that B? can be 
gauged away. 

In more detail we find that at first order 

i 

4 /2 = 0, £ - /2 = --e 1/2 , (5.9) 

with (1 + T r )ei/2 = 0. This data allows us to continue to three more orders and we 
find 



e+ = -§^i/2, e_ 1/2 = , (5.10) 



and 



and 



£± 3/2 = -igei/2, £_ 3 /2 = , (5.11) 

( i 2 i 2 c\ 

At the next order we find a new gauge parameter that has positive chirality. The 
only information that we require is that it can be used to gauge away B r , so we shall 
not write down the next order of the expansion explicitly. 

We can now use these results to determine the gauge transformations of and 
B & . For A & we find that 5A & = (p & — |r^)ei/ 2 . It turns out that B & is not 
the physical data to focus on. The asymptotic expansion of the gravitino equations 
can be specified by the data A& and B& or alternatively A& combined with a linear 
combination of B& and A&. We will fix this ambiguity, uniquely, by demanding that 
the latter data is gamma traceless, conserved and gauge- invariant. In terms of the 
symmetry algebra of the asymptotic AdS space, this decomposes into independent, 
irreducible representations. In more detail we define 



where 



and 



B & = B & + \iZ0A & + M^Ap , (5.13) 



M & p = P & %,P a ,, pJ = 5i-\T & T^ (5.14) 



3i£ 3i£ 3i£, 

3i - 9i 

;C^ & T^p0 ) + —-C^p a p &P 0. (5.15) 
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We then find that satisfies 



1 - T f )B & = 0, T & B& = 0, p & B& = 0, 6B & = 



(5.16) 



as advertised. It is convenient to also record here the constraints satisfied by the A& 
data and their gauge transformations 

(1 + r ) A & = 0, T & A & = 0, 5A & = (p & - \T & f) e 1/2 . (5.17) 
5.2 Extracting the Green's function 



The asymptotic expansion is specified by A& and B& satisfying (5.16 ), (5. 17) with, in 
particular, 



= r A ~ d A & + ■■■ + r~ a B & + ... 



(5.1* 



where A = 5/2 and d = 3. Using the expression for the D — 4 gamma-matrices in 
terms of the d = 3 gamma matrices presented in appendix A, we can write 

A & = ^ 2 (0, o a ), B & EE £ 9 / 2 (6 a , 0), £ 1/2 EE ^(O, e), (5.19) 

and the factors of £ have been added to give the canonical d = 3 dual field theory 



dimensions (we deduce from (3.1) that ip(r,p) has mass dimension —5/2 and recall 



the comments after (2.5)). The conditions (5. 16), (5. 17) can then be written in terms 
of the d = 3 spinors a a , b a and e as follows: 



7 a a a = 0, 5a a = (5f - \l a l P ) Ppt 
-1 a b a = 0, p a b a = 0, 56 a = . 



(5.20) 



Since the supercurrent is an operator in the dual SCFT with scaling dimension 
A = 5/2, the source should be fixed by the r _1//2 expansion data and the expectation 
value by the r~ 5//2 expansion data. After comparing (5.20) with (2.2), ( |2.4 ), the above 
expansion allows us to identify a a and b a as the source and the (gauge-invariant) 
expectation value of the supercurrent in the dual d = 3 field theory, respectively. 



Furthermore, (2.1) allows us to write 



G a fiO? = it 



(5.21) 



where ea\ i — 1,2, are the basis of 3d vector spinors that we introduced in section 
5} Recall that in section [2] we argued that the 0(2) symmetry of the system can be 
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used to show that t\2 — t 2 \ = and that t 22 (u, k) = t n (u, —k). Using these results 
we can obtain the convenient formula 



' 11 



ie { a ] b a 



(5.22) 



We can actually derive these results directly. In doing so we will obtain another 
formula for tu that will make contact with the gravitino equations that we derived 
in section 4 that involved the components and ip?. Recall that we decomposed 



the gravitino into 77 and p sectors (see (4.3), (4.5) and (4.7)) and that the resulting 



equations in the two sectors do not mix (in the 77 sector the equations for are 
given explicitly in (4.15) and in the p sector the equations for the can be obtained 



from (4.15) after taking k — >■ —k). To proceed we use the constraints (5.20) to write 

bf, =-- I 7 



Q 1 
~k at+ k P ' <1 ' 



dy 



7 



ty 



■<7/ 



1 



-T 9 * 1 at + ^-p-a. 



(5.23) 



We can also use (5.3), which says that asymptotically we have 

2i 



p ■ a 



-a (1) 



(5.24) 



the superscript is a reminder that it appears at one lower order in the expansion in 



1/r. Then taking the t component of (5.21) we are led to 



b~ t = -2[t u n - t 21 m]n ( p 2 a t ~ - — (k - 2u)a 



2[t 2 2m - t 12 n]m p —{k + 2u)a~ 



(i) 



(5.25) 



Now, after imposing ingoing boundary conditions at the black-brane horizon (which 
we discuss in subsequent sections) we will have two solutions in the 77 sector and two 
solutions in the p sector which allow us to extract the four components t%j. Consider 
first the 77 sector solutions. Since Af and = ^.p ■ A both have negative chirality 
with respect to T r , they must be proportional to (r] + + 77") = (0, m) T in the basis 
we are using. Similarly, since has positive chirality with respect to T r , it must be 
proportional to (77+ — 77^) = (0, n) T . Thus, by considering the 77 sector solutions, and 
recalling rhn = —nm = 1, we deduce that t 2 i = and that 

mbt 



hi 



2n(-p 2 a t + Uk - 2Co)a { r ) ) 



(5.26) 
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By a similar chain of reasoning in the p sector we deduce that t 12 = and that 

nb t 



t 



22 



(5.27) 



2m(p 2 a t + j,(k + 2Cj)ar ] ) ' 

That £22(^5 k) = tn (a;, ~ k) follows from the fact that the equations for the in the 
p sector are exactly the same as those for the in the r\ sector after taking k — > —k 
and also that we have (p + + p _ ) = (0, n) T and (p + — p~) = — (m, 0) T (note the minus 
sign). 

To summarise, we can obtain t n , which contains all of the independent infor- 
mation in the retarded Green's function, by solving the rj sector gravitino equations 
for the given in (4.15), with ingoing boundary conditions, extracting the asymp- 
totic behaviour at the AdS boundary and then substitute into (5.26). It is worth 
emphasising the following point. The equations (4.15) immediately give the t and f 
components of the gravitino and hence a£ and a?. However, to obtain the correct bf, 
one also needs to consider the x and y components and ensure that the conditions 
(5.16) are satisfied. 



6 The Green's function for 

Having understood the boundary asymptotics of the N = 2 gravitino, which is in- 
dependent of whether the AdS-RN solution is extremal (T = 0) or non-extremal 
(T > 0) we must now address the near-horizon behaviour, which shows a dichotomy 
between the two cases due to the different singularity structure of the equations. In 
both cases, in order to determine the retarded correlation function we want to spec- 
ify ingoing boundary conditions. In this section we consider T > 0, analysing the 
near-horizon behaviour and then presenting the results of our numerical integrations 
for the Green's function t n , focussing on the spectral function A =Im t n . We will 
consider T = in subsequent sections. 



6.1 Near-horizon series 

At temperature T / the function / appearing in the AdS-RN black-brane metric 



(3.4) near the horizon has a simple zero: 

f(r) = 47rT(r - r+) + . . . (6.1) 
Thus, near the horizon the tortoise coordinate is given by 

u = ^_=iA, (6 . 2) 
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Furthermore, the simple zero in / means that the equations (4.15) have a regular sin- 
gular point at r = r + and we can apply a version of Frobenius' method to the system 
of four coupled ODEs. This is done most easily in the first-order formalism, where 
the independent solutions are determined as the leading-order eigenvalues and eigen- 
vectors of the matrix multiplying the non-derivative terms in the coupled equations. 



Concretely, we write the system of differential equations (4.15) near the horizon as 

-A(ou, k, r)u(u, k, r) — , (6.3) 



w d/2) 



dr r — 

where the function u = (u^ 3 ^ 2 \ w 1 ' 2 ', w* -1 / 2 ), w(~ 3 / 2 )) and A(tu, k, r) is a 4 x 4 matrix, 
one can show that we have a series solution in ascending powers of (r — T+) 1 / 2 , with 

oo oo 

A = An(w, k)(r - r + ) n ' 2 , u = ^ u n (u, k)(r - r + ) n / 2 . (6.4) 

n=0 n=0 

The leading-order result for the matrix A(u, k) reads 

Ao = dia H~4-4^T'4 + 4^T'4-4^T'-4 + 4^Tj ' (6 ' 5) 
Being diagonal, we immediately obtain the indicial roots and the associated eigen- 
vectors and hence the four different solutions at the horizon. 

Since ingoing solutions should behave like e~* tJ (* +r *) near the horizon, we deduce 
that we should impose the horizon boundary conditions 

= & 4 V2) = 4~ 3/2) = • ( 6 - 6 ) 

•+ 

The two independent ingoing solutions then have leading-order behaviour 

u w = (r-r + )-H? (4 s) + 0((r-r + ) 1 / 2 )) , s = |, -| (6.7) 

where we have introduced the often-employed dimensionless measure of frequency at 
finite temperature: 

(s) 

and the u are two arbitrary complex coefficients. Note that the near-horizon expan- 
sion is a series in ascending powers of (r — r + ) 1 ^ 2 and that these two free parameters 
determine the entire series expansion of u (all four components are non-vanishing) to 
arbitrarily high order in principle^} Gauge fixing restricts this to a single physical pa- 
rameter in this sector as we shall see now. In particular, this parameter arising in the 
r) sector of the gravitino and the corresponding parameter in the p sector correspond 
to the two physical polarisation states of the gravitino. 



5 And to very high order in practise, using computer algebra. In our numerical procedures (see 
appendix C) we expand to order (r — r + ) 15 / 2 to set the boundary conditions near the horizon. 
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6.2 Near-horizon SUSY 



We now consider the residual supersymmetry transformations (4.16), (4. 17). Writing 
w = (m/ 1 / 2 ), w/ -1 / 2 )) the equations (4.17) have an expansion 

dw(u, k, r) 1 



dr 



r — r i 



-B(oj, k, r)w(oj, k,r) = 



(6.9) 



where, as before 



(6.10) 



n=0 



n=0 



This time find the leading-order matrix 



1 iu 1 iu 
B = diag 7 - ^ , 7 -+ 



(6.11) 



4 4ttT ' 4 4ttT, 
and thus, the ingoing gauge transformation satisfies the horizon boundary condition 



w 



(-1/2) 



& w { 1/2) = . 



(6.12) 



We can then develop the most general ingoing gauge transformation parameters as a 
series with leading form 



w 



(1/2) 



(6.13) 



with the arbitrary complex number w^ 2 ^ determining the entire series for w(oj, k, r). 
Plugging this into the SUSY variations (4.16) finally gives the desired expressions for 
the leading ingoing pieces 



5u {s) 



r — r_i_ 2 2 



(/«(»,) + ."), 8 



2 ' 2 



(6.14) 



where 

/ (3/2) M 
/ ( - 1/2) M 



—w, 



w, 



(1/2), 


( i/ 2 ) 1 2kH 2 -2k(e 2 jj + ir + ) +£ 2 jj 2 + 4ir + (jj + Lo) +r 



/ 7rT(| + ztt)) 



.2 «-2 



4r 2 VvfT 



«tt> 



(6.15) 



From this calculation we can see that the two parameters defining the near-horizon 
ingoing solutions transform under the residual gauge transformations as 



^ S) = /W(m) 



(6.16) 
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Figure 1: Representation of spectral density A(u,k) as a function of frequency and 
momentum at T/fi = 0.04. Panel a) shows a surface plot in the (ou, k) plane with u- 
slices taken in intervals of A(k£) = 0.5. Panel b) shows an alternative representation 
where a selection of such u slices for positive ki e (0.1, 1.1) is shown superimposed on 
one plot, with larger values of kl in darker shades of grey. Furthermore the scale of 
the peak at kt = 0.1 is shown separately in an inset. The gap between the two regions 
of high spectral density is clearly visible in the first panel (for more information see 
Fig. 0). 



which can be used to set one of the to zero. Alternatively one can focus on a 
single gauge-invariant combination of initial conditions 

(-1/2) / ( ~ 1/2) M (3/2) 

< ~ /( 3/ 2)(m) u o • (6-17) 

In practise, for our numerics, we always choose a specific gauge, but have checked 
extensively that the resulting expressions extracted for the correlator are manifestly 
independent of gauge. This is not only reassuring but also gives a strong check on 
our numerical algorithms. 

6.3 Numerical results at T > 

Using the numerical procedure outlined in appendix [D] we can now solve the gravitino 
equations in the r] sector and then extract tu(u, k) from ( |5.26 ). We are especially 



interested in the spectral function A(uj, k) = lmtu(u, k) which encodes information 
about the density of states at finite temperature. A surface plot of the spectral 
density is shown in Fig. [T] together with another representation of the type used in 
the companion paper [l] . A collection of other ARPES-likeU>lots showing the spectral 

6 ARPES, or angle resolved photon emission spectroscopy, is an experimental technique in which 
the charged (surface) excitations of a material are measured and represented as a density in the 
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density A(u, k) as a function of u and k are shown in Fig. [2j Regions of high spectral 
density are displayed in 'warmer' colours. 

The most striking feature on the (u, k) plane is the region of high spectral weight 
around k = 0,u> = —fi. This long-wavelength mode is associated with the phonino, a 
collective fermionic excitation of any supersymmetric medium at finite temperature 
and density, that we discussed in section 2.2. The pole is at k — but frequency 
oj = —fi due to the presence of a finite chemical potential, which arises as a weak 
gauging of the i?-symmetry of the boundary theory. 

There is a second region of high spectral weight, albeit suppressed with respect 
to the phonino excitation, located at positive k. For all k there is a depletion of the 
density of states near the origin, with A(ou, k) getting closer to zero around u = 
as the temperature is lowered. We will see in later sections that at strictly zero 
temperature, the spectral weight actually vanishes at the origin and that there is a 
power-law soft gap in the spectrum. 

6.3.1 Supersound propagation and diffusion 

Returning once more to the phonino, that is to small k and small u, we assume that 
the retarded propagator takes the form dictated by the presence of a simple pole 
located at u = k = 0: 

u> — il[k, /i, T) + iY[k, fi, T) 

with Z, Q, T real and Z and F positive at k — 0. A straightforward analysis of 
A(u, k) = Imtii(<D, k) shows that the spectral peak at u = Q has height Z/Y and 
that the full width at half maximum is Ad) = 2r. 

Given the above analysis, our numerical results allow us to fit a dispersion relation, 
encoding the location of the maximum of the peak as well as its broadening (its full 
width at half maximum), of the form 

rb = v s q-iZ(T,n)q 2 + ■■■ (6.19) 

where q = k/2irT and to the given accuracy we fincQthe propagation speed 

w s ~ 0.50, (6.20) 

(w, k) plane. 

Note that there are two sources of error here. Firstly numerical, but secondly also the error due 
to the fact that the dispersion itself is only true in the sense of a power series for small values of to 
and q. 



23 




-1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 
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Figure 2: Spectral density A(u,k) for high and low temperatures (in units of the 
chemical potential): T/fi = 0.37 for panel a) and T/^i = 0.04 for panels b) & c). In 
panels a) and b) values of low spectral density are in shades of blue, whereas regions of 
high density are in shades of red and the largest peak is cut-off and the cut-out region 
shown in light pink. At finite temperature the most prominent feature is the large 
spectral density region due to the phonino, or supersound pole, in the left bottom part 
of the (u , k) plane. There is a second region of larger spectral density in the positive- 
frequency half plane and as the temperature is lowered a region of suppressed spectral 
weight of size O(n) opens up between these two regions, that is approximately power- 
law. However, as illustrated in c), for momenta kl G (1.2,2.1), thermal excitations 
cause A{oj <C /i, k) to be nonzero at any T^O (see also zoomed region around origin 
in inset). However, A(u),k) can be very small near the origin. 

and we shall return to the diffusion coefficient presently, after giving an analytic 
argument for the constant speed of propagation v s . Indeed our fitted value for v s is 
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consistent with the value v s = 1/2 that is expected for conformally invariant field 
theories 65 . In fact we can obtain some further insight as follows. The normalisable 
phonino mode should have A a = in the asymptotic expansion of the gravitino. 
In an expansion in small values of u and k we find that at leading order there is a 
solution to the gravitino equations with = C & &T"e, where e is a constant spinor 
and the diagonal matrix C was given in (5.7). This can be viewed as a large gauge 
transformation and satisfies T a Ba = 0. The condition p a Ba = then becomes the 
Dirac equation (2u)r* + kT x )e = with phase velocity v s = 1/2. 

Now returning to the diffusion coefficient, we note that in the absence of a chemical 
potential, simple dimensional analysis would dictate that £ = const. One then 
defines the diffusion constant D s = 27rT£, governing the broadening of the peak 
in u, k space. Since in the present context we have a non-zero chemical potential, 
dimensional analysis is no longer sufficient to ensure that D s /(2ttT) is a constant, 
but rather it allows for it to be a non-trivial function of the dimensionless ratio of 
temperature and chemical potential, Ds ^^ ■ 

Our extensive numerical data allows us to determine D, for a number of values 



of T I fx and we show the results in Fig |3j Note that it approaches a constant in the 
appropriate limit T 3> /i, where the zero chemical potential result should apply. It 
would be interesting to carry out a more exhaustive analysis, which is however, beyond 
the scope of the present work. Finally, it is worthwhile to point out that the same 
dimensional analysis argument would in principle allow for the propagation speed 
v s to be a non-trivial function of T/li, but the additional requirement of conformal 
symmetry, which is at heart of the brief analytical argument we gave above, ensures 
that it is actually a constant, in complete agreement with our numerical results. 



We should compare the results of this section to the work of 67 which used an 
analytical expansion at small frequency and momentum, but zero chemical potential, 
to extract the analogous information in five bulk dimensions where it was found that 
v s = 1/3 and £ = A^/2/9. Note, that due to the absence of finite /i the diffusion 
coefficient is a constant, as expected from dimensional analysis. 
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Figure 3: Supersound diffusion rate. For high values ofT/fi the diffusion coefficient 
approaches a constant, which is expected since in this limit it should approach the 
/i = form and by dimensional analysis must be a constant. 
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7 Extremal T = case 
7.1 Near horizon series 

At zero temperature the horizon at r = r + is an irregular]^] singular point of the 
gravitino equations, arising from the confluence of the two regular singular points at 
r±. However, we are still able to develop a systematic expansion by factoring out the 
leading essentially singular behaviour. The remaining equations then, again, allow a 
Frobenius series expansion, as we describe now. The extremal limit of the function 



/ appearing in the metric (3.7) has a double zero at the horizon 



f~Lf 2) (r-r + y + ... (7.1) 

which gives rise to a confluent singularity in the gravitino equation. The tortoise 
coordinate near the horizon is given by 

n = @- • (7.2) 



8 At zero frequency this irregular behaviour is absent and we can develop an ordinary Frobenius 



series. We can sidestep this subtlety for now, but will return to it in section 9.3 where it will allow 
us a valuable analytical insight into the spectral function at zero frequency. 



26 



-1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 

a) 




Figure 4: a) Spectral density A(u, k) at T = 0. As before, the largest peak associated 
with the phonino pole, is cut-off and the cut-out region shown in pink to distinguish 
it from the white areas of the plots. Values of low spectral density are in shades of 
blue, whereas regions of high density are in shades of red. b) The spectral density 
at discrete values of ki £ (0.1,1.1), with increasing values of momentum in darker 
shades of grey, c) For higher values of ki £ (1.2,2.1), as is already apparent from 
panel a), the bump at positive frequency becomes the dominant feature of the spectrum. 
The inset highlights the emergence of a power-law gap at zero frequency by fitting the 
num erics close to the origin to the analytic result oc u 2uk , where as shown in section 



9.2 



v k - \l Y2 + 2j? 



After defining the quantity 



A(w) * + ^ = J_ + S «, (7.3) 
2^3 3r+ 2^3 3^3 H 
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we find that there are four independent solutions, specified by the four complex 



coefficients with the following leading behaviour 



r 2 

= e"^V (r - r+ ) s+iA ^ + C(r - r+)) , s = §, -§ . (7.4) 

Note that the leading essentially singular behaviour is factored out. Also note that, in 
contrast to the finite-temperature solution, the series solutions proceed in ascending 
powers of (r — r+) (rather than its square root as we saw for T > 0). Since the first 
two solutions behave in terms of the tortoise coordinate as e~ tur * , we identify them 
as ingoing at the horizon, whereas the latter two are outgoing. Hence, also at zero 



temperature we impose the boundary conditions u^ 2 ^ = Uq 3 ^ = 0, as in (6.6), 
leaving two solutions. 

The parameter of the residual gauge transformations can be similarly factored, 
allowing us two write down two independent solutions, specified by the complex 
parameters w^ 1 ^ 2 ^, with leading behaviour 

.,2 

w (i/2) = e ^^(r-r + )-^- iA ^(w^ /2) + 0(r-r+)), 



2 



-iL 

w (-m = e ^T( r - r+ )-|+i A H(4- 1 /2) + ( r _ r+) )_ ( 7 _ 5 ) 
We thus impose Wq l ^ 2 ^ = as the ingoing horizon boundary condition. Substituting 



the resulting expansion into the gauge variation equations (4.14) gives us 

§u (s) = e -^ g (s)^ w (l/2)( r _ r+ )-»- iA( W ) + Q( r _ r+ y ? ^ 7 _g) 

for s = |, — | , where the g^ s \uj) are helicity-dependent constants, the precise form 
of which will not be needed. Again, we see that by virtue of the residual gauge 

(3/2) (— 1/2) 

transformation we can set either of u = or u = 0. As at finite temperature 
there is a gauge-invariant combination that follows from the action of the gauge 
transformation on the horizon parameters, but since its exact form, the exact analog 



of (6.17), is not needed and not very illuminating we do not present it here. Thus, as 
expected, gauge fixing restricts us to a single physical parameter in the rj sector and, 
similarly, a single parameter in the p sector and these correspond to the two physical 
polarisation states of the gravitino. 

7.2 Numerical results at T = 

Using the numerical procedure outlined in appendix [D] we can now solve the gravitino 



equations and extract t\\{u, k) using (5.26). Recall, once more, that knowing tu(uj, k) 
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for positive and negative values of the momentum also gives us £22 k). A plot of 
the spectral density A(co, k)=lm tn(u, k) is shown in Fig. |4j 

We find a structure that is very similar to the finite temperature results, as one 
might have expected by continuity. The phonino pole is still responsible for the 
dominant feature on the (u, k) plane, although of course its interpretation in terms of 
hydrodynamics is rather subtle at zero temperature. A discussion of this issue would 
go beyond the scope of this paper. Again, there is a second region of high spectral 
weight at larger positive values of momentum, separated from the phonino peak by 
a power-law gap. Interestingly, at zero temperature, the spectral density actually 
vanishes at to = exhibiting a characteristic power-law oc tu 2uk , with = y^~+ ^2, 
for small values of frequency, as shown clearly in Fig. [4)3. As the exact expression 
for z/fc suggests, we can in fact derive this power law analytically as a low-frequency 
expansion of the retarded correlator. This expansion makes heavy use of the exact 
solution of the gravitino equations in the near-horizon geometry of the extremal limit, 
AdS2 x IR 2 , of the AdSRN black-brane, which we will turn to now. 



8 Analytic solution in the AdS2 x JR 2 background 

In this section we pause to solve the gravitino equations exactly in the AdS2 x M. 2 
background. In the next section we will use these results to analyse the non-analytic 
behaviour of the spectral function in the extremal AdS-RN black-brane background 
for small u. 

The AdS2 X K 2 background is given by 

ds 2 = \-dt 2 + dz 2 } + r -±dx l dx\ A = -®-dt , (8.1) 
z 2 1 J £ 2 y/2z 

and we note that the boundary of the AdS2 is at z = and the Poincare horizon is 

at z = 00, where z = L 2 2 ^/(r — r + ). In the frame 

e ± = [_ dz ± dt j ? ( 8 .2) 
we find tt± = =Fl/(4L( 2 )), A± = ±1/(2^2), / = -1/(2V2L (2) ), and hence the 



gravitino equations (in the rj sector) (4.8) read 



Al3« (3/2) - K+ k u^ = 0, 

A+uW*> - K+ h uW> = 0, 

AtJ- 1/2) -K-u^ = 0, 

A+. ±u (-W _ K - u (-im = 0, (8.3) 
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where we have defined the differential operators 



.4 



± ._ 



zd z + n ± z 



UJZ + 



g£(2) 



and also 



K i = £(2)Mfc ± 



71' 



ji k = m + 



iki 
r+ 



U) 



$.5) 



Notice that in (8.3) the Pauli coupling essentially contributes a shift of the mo- 
ntum (this was also observed in bottom- up models for spin 1/2 fe 
Similarly, the residual supersymmetry transformations (4.14) read 

5u^ 



mentum (this was also observed in bottom-up models for spin 1/2 fermions in 38 43]). 

asymmetry 

— Atw^ , 

2 

A\w^-L {2) 

2 

A~_M l M - L {2) ( m + 

a-M- 1/2) 



5uW 
5u^ 



2L( 2 ) 2 
1 

2L {2) 
1 

2L {2) 
1 



m — 



\/2L( 2 ) 



(1/2) 



(-1/2) 



(2) 



2L 



(2) 



.6) 



where, from (|4.12|), the residual gauge transformation parameters satisfy 
A-_ lW ^-L (2) (2m+— V" 1/2) = 0, 



a + M- 1/2) 



2m — 



r+ J 
ik£\ 



- w 



(1/2) 







(8.7) 



From (8.3) we can derive the following second-order equations for the component 
functions 



-K~ Kn 

V I -I ^ fc « 



u 



M (3/2) 

(-1/2) 
(-3/2) 



0. 
0. 

o. 

0. 



These equations can be solved exactly. Before doing so, let us first analyse the 
asymptotic behaviour near the AdS2 boundary z — > 0: 



u (- 1/2 ) )U (-3/2) 



Az*- U+ (1 + 0{z)) + Bz* +V+ (1 + 0(z)) , 
Az*~ v - (1 + O(^)) + Bz^ +U - (1 + , 



■9) 
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where v± 



1 + L\ 2) 



|A 2 



1/2 



It is interesting to 



|/i fc | 2 - \q 2 ) ±iV2L {2) (q + m\) 
observe that the exponents v± always have imaginary parts, unless q = —mX, which 
is precisely the case for the supersymmetric N = 2 theory we are interested in. Indeed 
in this case the two exponents in fact coincide: v± = Vk with 




7 k 2 
12 + 2jj? 



(8.10) 



The expansion (8.9) implies that the scaling dimension of the operators, labelled by 
the spatial momentum k, in the dual one-dimensional CFT, is given by 1/2 + v^. 

We now see that there is no possibility for the log-periodic behaviour that was 
seen in the bottom-up models for spin-|in pj], since this feature only arises when 
the IR conformal dimensions are not real. The reality of is due to the charge to 
mass ratio in minimal N = 2 supergravity and this can be viewed as being uniquely 
fixed by the structure of the embedding in string/M-theory. From a physical point 
of view the same ratio also means that pair production of charged Fermions is not 
energetically favourable in the constant electric field in the AdS2 region jH]. This 
suggests that the system does not exhibit an instability towards forming an electron 



star [18j|35j|53]. Note that in (8.9) we have labeled the source behaviour with the 
letter 'A' and the expectation value with l B\ 



We can in fact write down the general solutions of the equations (8.8) the asymp- 
totics of which we just explored. The second-order equations all follow the template 



z l di + 



2ia(s)uz + {wzY 



u. 



0. 



where 



a(s) = s + 



V2 



(8.11) 



(8.12) 



This is Whittaker's equation and there are two linearly independent solutions which, 
for each vS s \ we can group into 



awW a ,u k {-2iuz) + aovrW- ajVk (2iuz) 



(8.13) 



where W av . is Whittaker's confluent hypergeometric function. The first term in 



(8.13) corresponds to a purely ingoing solution, whereas the second term gives the 
purely outgoing one^j This can be seen from the asymptotic behaviour for large x 

{! + ■■■). (8.14) 



W a ,v k (ix) 



} Notc that in the extreme AdS-RN black-brane, near the horizon at r 



the tortoise coordinate 



is given by 



r— r+ 



-z. Thus the ingoing boundary condition is e 



-icj(t+r w ) 



-iu{t — z) 
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It is now easy to write down the ingoing and outgoing solutions. As we are most 
interested in the retarded Green's function we will only show the ingoing solutions, 
for which = a^W a ^ s ^ Vk (—2icoz). The first-order equations (8.3) will now relate 
the constant coefficients a^ s \ so that only two free complex parameters remain. Using 
the identities 

= -(\t-a)W aiVk {t)- [i/g- (a- |) 2 ] W a - 1>Vh {t), (8.15) 
we find a/ 1 / 2 ) = — a^'^K^ , a _ ( 3 / 2 ) = — a^^KZ , and so we can now write 



M (3/2) 



0(3/2),*/* 



-2iujz) 



-a^K; k W a{l/2) , Uk (-2iuz) , 



-2iujz) . 



(8.16) 



Analogous expressions hold in the p-sector. After substitution in (4.7) and then in 
(4.3) we obtain the ingoing behaviour of the gravitino in the near-horizon AdS 2 x M. 2 
region. 

Finally we consider the residual supersymmetry transformations. We find that 



we can again solve (8.7) using ingoing Whittaker-type functions 

w (1/2) = 2d 1 W a(1/2hUk {-2iuz) , 
w t~W = 2d 2 W a{ _ l/2>k (-2iuz) , 

where q are constants satisfying 

'ikl 



d 2 = L 



(2) 



2m d\ . 



(8.17) 



(8.18) 



After substituting into (8.6), and also using the supersymmetry values for q, m, A, we 
obtain 



W, 



'(2) 



a(3/2),u k 



-2iuz) 



L <i) 

do 



K+ k W a(1/2) , Uk (-2iuz) 



K n k w <*(-i/2),v k {-2iuz) 



/ —2iuz) . 



.19) 



32 



Comparing with (8.16) we see that the residual supersymmetry transformations can 
be used to set either w^ 3 / 2 ) = i^ 1 / 2 ) = or alternatively u^ 1 ^ 2 ^ = w^ -3 / 2 ) = 0. 

By expanding the solutions (8.16) near the boundary of the AdS2 x IR 2 region 
we can extract the retarded Green's function of the dual one-dimensional SCFT. 
Specifically, introducing the constant vector-spinors 



si 



S2 



s 3 



s 4 



® rf - T 



a(3/2) + v k 
Kt 



i _ a (3/ 2 ) _ „ fc 



e ® ?7 — y 



-a(-l/2) + i/ fc 



e <g) 77 



e Cg) 



a(-l/2)-i/ fc 



(8.20) 



we find that after substituting (8.16) into (4.7), in the limit that z — > 0, we have 



X 



+ a 



(-1/2) 



si + £ R (w,a(3/2),z/ fc ) 



'(2) 



s 2 



-"k 



s 3 + G R (u,a(-l/2),v k ) 



«4 



(8.21} 



and the bars on a^ 3 ^ 2 \ aS 1 ' 2 ' indicate that we have absorbed some overall factors. 
In these expressions 

T{-2v k ) V{\ - a(s) + u k ) 



(2uL {2 )) 2Uk 



T(2v k ) Ta-a(s) 



v k , 



(8.22) 



define, up to normalisation, the retarded Green functions of the emergent IR CFT for 
each helicity state, that is each one-dimensional representation of Spin(l, 1) labeled 
by s. We show more details on this computation, as well as its cousin with advanced 
boundary conditions, in appendix [E| Note that the residual supersymmetry allows us 
to choose either the gauge a*- 3 / 2 -* = or a^~ 1//2 - ) = if we so wish. 



9 Matched asymptotic expansion 

In this section we derive the leading-order result of a small-frequency expansion of 
the retarded correlation function at zero temperature and finite chemical potential. 
This leads to an explicit expression for the scaling properties of the spectral function 
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at small frequency. The analysis is very similar to the one first used for spinors in [5] , 
but we face additional complications stemming from gauge invariance as well as the 
multi-component nature of the gravitino. Our treatment follows closely that of [H], 
but our emphasis differs in some respects. 

We start by discussing carefully the near-horizon limit in which the analytic so- 
lution in AdSi x IR 2 is applicable, and then how this can be matched to the leading 
small-frequency solution at large r. 



9.1 Scaling limits and the small u expansion 

When T = the only physical parameter in the bulk solution is the chemical potential, 
given by \ii = y/3r+/£. Naively the solutions in the AdS 2 x M 2 background given 
in section M are applicable close to the horizon, as measured relative to the chemical 
potential, that is for (r — r + )/r + <C 1. This is the regime in which the higher-order 
terms in / and in (3.7) and (3.5) can be dropped. However, this limit requires 
some care since the term u/f in (4.15) diverges faster at small (r — r+) than the 



other terms. To get the precise AdS 2 x 
small-frequency limit where cj//x <C 1 • 



equations (8.3) one must also take a 



To see this explicitly, let us define the new coordinates 

e T 2 



r — r i 



(9.1) 



so that taking the limit e — > 0, keeping £ and r fixed, of the metric and gauge field 
(3.4) gives the exact AdS% x M 2 near-horizon geometry. In other words 



ds 2 



T 2 

C 2 



(-rfr 2 + d( 2 ) + r ^dx i dx i + 0(e) , A 



(2) 



2C 



dr + 0(e) . (9.2) 



We now perform an expansion of the gravitino equations (4.15) in the parameter e 
to find that (at zero T) we get the set of equations 



e 

^(1/2) + 

e 



3 iqL m 



ti 



(3/2) 



1 iqL {2) 
2 + V2 



(3/2) 



e 
e 



C5 C + 



V2 J 



■u 



(-1/2) 



(-3/2) 



[J2 ) + 0(e) = 

-K£ k uW + 0(e) = 

K- h + 0(e) = 

W (-V 2 ) + 0(e) = 



(9.3) 
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The frequency u appearing here is the Fourier conjugate to the time coordinate t, 
but in order to obtain a regular near- horizon limit, we were forced to scale et — r and 
thus the relevant frequency for the near-horizon region is instead u T = u/e. With 
this identification the leading-order behaviour is actually O(e ). The equations at 



leading order then precisely reproduce the AdS-z x ^ 2 problem (8.3) we solved above, 
now written in terms of ( and r (or rather its Fourier conjugate u T ). 

Hence to get the well-defined set of equations for finite u T and (, we want to take 
the limit 

u£ ->■ e -> , with u T , ( finite . (9.4) 
Recall that the corrections to the metric function / and potential (f> are small provided 



(r — r + )/r + <C 1. Given the coordinate redefinition (9.1), the limit is thus a good 
approximation provided 

< r ~ r+ < e, u£ < e, (9.5) 

with the size of e fi£, set by the chemical potential. 

In the language of matched asymptotic expansions [3], this AdS2 X limit de- 
scribes the inner region for a perturbative expansion of the differential equations in 
small uj. It is subtle because we are considering a range of r where the u/f term is 



diverging faster than the other terms in (4.15 ), and so cannot simply be dropped even 



though uj is small. The outer region is considerably simpler. It is given by the range 
of r where the u/f term can be safely ignored relative to the other terms. Since 1/f 
diverges as t/(r — r + ) relative to the other terms, the outer region is defined by 

e'<^+<oo, (9.6) 

with e' ^> u£. In each region we can expand the solution as a power series in u£ 

^ = ^ 0) (r) + M)4 1) (r) + -- - , 
^o = ^\r) + {u£)^\r) + ... , 

where the subscripts "I" and "O" denote inner and outer regions. 

In this case the matched expansion is very straightforward. If we take u£ small 
enough we can simply make the two regions overlap and match the solutions. Specif- 
ically, by taking a double-scaling limit 

u£ < e < e < fi£, (9.8) 

both inner and outer solutions will be valid in the region e' < (r — r + )/£ < e. Let us 
now turn to the resulting leading-order matched solution. 
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9.2 Leading-order solution 



By comparing (9.3) with (8.3) and (8.16) it is clear that the leading-order (in uj 



ingoing solution in the inner region takes the form 



(z; uj, k) oc W a ( 8 ) )Vk (-2iuz) 



(9.9) 



with the exact coefficients of the individual components as given in (8.16) and where 
we are using the coordinate z = L%J (r — r+) = £ 2 /6(r — r + ). In the matching region 
we have 

w£<e', (9.10) 



e < 



z < e 



and hence wz < 1. Thus we can use the asymptotic, near-boundary expansions 
(8.20) from the AdS^ analysis of section [8] Explicitly we have 



Xl ~ a {3/2) 



T 2 
r — Ta 



i-Vh 



si + G R (u),a(\),v k ) 



T 2 
r — r i 



S-2 



(-1/2) 



+ a 



T 2 

r — r a 



s 3 + G R (uj,a(-±),v k ) 



T 2 
r — r a 



(9.11) 



Note that Gr{uj, ot(s), U}.) are non-analytic in u. In addition, we could have fixed the 
gauge, e.g. by setting a^ 3//2 ^ = 0, but we prefer, for the time being, not to do so. 

The leading-order outer solution is obtained by solving, at least formally, Eqs. 
(4.15) at cu = 0. As was already observed in |3J for spin 1/2 fermions, at zero tem- 



perature and strictly zero frequency the essentially singular nature of the horizon is 
absent (the problematic terms are proportional to uj) and we can develop an ordinary 
Frobenius series. The leading terms in this expansion give the approximate outer 
solution in the overlap region. Prior to gauge fixing there are four independent solu- 
tions and these can be fixed by their asymptotic behaviour near the horizon. For this 
purpose we can define X(i) with i — 1,2,3, 4 such that for e' < (r — r + )/£ < e \it 
we have 



X(i),(3)(r, k) 



T 2 

r — r i 



s l,3; 



X(2),(4)(r, k) 



T 2 

r — r i 



S2A- (9.12) 



As expected, we see that this behaviour is identical to the small-argument asymptotics 
of the inner region. Matching to the inner solution fixes the coefficient of each xu) 
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component and we have 

Xo(r,u,k) = a (3/2) X(i)(r, k) + Gr(co, at(§), Vk)X(2)(r, k) 



+ a 



(-1/2) 



X(3)(f", fe) + Gr{u, Qf(-i), ^fc)X(4)(^, fc) + CM). 



(9.13) 



More generally, recalling the asymptotic expansion of the Whittaker functions (E.l ) 
we can extract the non-analytic w-dependence of the outer solution by writing 



Xo(r,u,k) = x v {1) (r } uj,k) + G(uj,is k )xl 2) (r,uj,k) . 
where x%\ are analytic in uj, with 



(9.14) 



X v {1) (r, u>, k) = a (3/2) X(i)(r, k) + a^ 1/2) X(3)(r, k) + 0(u£), 

Y V ( L.\ _ 3 (3/2)^ / LA , -(-1/2) K ~ 27l) K ~ 2?3 ~ X ) 



X (4 )(r, k)+0(oj£), 
(9.15) 



and we have defined 



Q(u, v k ) = Gr(u, a(§), i/jfc), 



-»7rr/ fc 



,r(-2^) r (- 1 - 2 73 + ^) 



r( 2 ^) r(-i-^-u k : 



(2uL {2 )) 



1v k 



(9.16) 



Note that our choice in factoring out G(oo, is simply one possible option corre- 
sponding to Gr(cu, ot(s), Vk) with s = 3/2. Any other choice of s, related to this by 
the familiar recursion identities of the Euler gamma function, would be equally valid. 
Furthermore all choices lead to the same non-analytic u dependence oc (uL(2)) 2v of 
the end result. 

We now have all the ingredients to write down the answer for tn(u, k). To leading 
order in the expansion we have 



tu{u,k) 



ie {2) 



(9.17) 



and of course the analogous expression for £22 (w, k). In these expressions the quanti- 



ties bJF and o^'" are the various near-boundary data (c.f. (5.22)) of the asymptotic 



"(<) 



expansion of the two solutions (9.15) at n-th order in the ul perturbative expansion. 



We have no analytic expressions for these coefficients as functions of k, but they can 
be determined numerically from the outer region equations. 
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9.3 Scaling behaviour of spectral function at zero frequency 

The immediate utility of the exact expressions obtained in the last subsection is 
that we can use them to give a simple analytic argument for the scaling properties 
of the spectral function near the origin. We shall establish the scaling property by 
combining the action of a discrete symmetry with the specific form of the horizon 
boundary condition at zero frequency and zero temperature. 

The discrete symmetry of interest is the three-dimensional inversion (t, x, y) — > 
(— t, —x, —y). This can be viewed as the combination of time-reversal symmetry T 
and the rotation R : (t, x, y) — > (t, —x, —y) already discussed in section [2j Here we 
will focus on the action of these symmetries on the bulk gravitino field. 

The action of time-reversal is 

T : ^(t, x, r) ^ ^\t, x, r) = -B T P/VvH, x, r)* , (9.18) 

where V = diag(l, —1, —1, —1) is the parity inversion matrix in 4d and Bt induces 
the transformations T* i-> B^^Bt = — T* while leaving all other gamma matrices 
untouched. The matrix Bt is uniquely defined up to an arbitrary phase, which here 
we fix so that Bt = T rx y . The complex conjugation is included in line with the usual 
action for solutions of complex wave equations which include a coupling to a U(l) 
field. It ensures that the bulk gravitino equation is time-reversal covariant given the 
usual map of the gauge potential *4 M (£, x, r) i-> V^ u A u (—t,x.,r) under T. From the 
boundary field theory perspective, where ip ^ is related to the expectation value of S a , 
it is a consequence of the anti-linearity of the time-reversal operator. 

Recall that the black-hole background is static and purely electrically charged 
(only Ai non-zero) and hence is invariant under time-reversal . Hence T should be 
a symmetry of the bulk gravitino equations of motion. To see this, note that the 



corresponding action on x^(t,r, k) defined in (4.3) is 



T Xi (t,rM)\ /-r«^,(- t , r ,-k). 
1 Yf(i,r,k / \ r"» x ,(-t,r,-k)« 1 



Expanding in the basis (4.7), this leads to 



u 



is) 



(t,r,k) i — y v(~ s \-t,r, -k)" 



T : { v ' ' ' v ' ' ' . (9.20) 

v^(t,r,k) h-> -u { - s \-t,r, -k)* 

Given time-dependence of the form e~ luJt , the frequency u is unchanged under this 



map, and we see that (9.20) is manifestly a symmetry of the bulk equations (4.15). 
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The action of the rotation R on the gravitino can similarly be written as 



(9.21) 



Here 1Z is the rotation matrix TZ = diag(l, —1, —1, 1) and the matrix T x ^ effects the 
appropriate rotation on the spinor indices. Acting on the components we have 



R : 



u {s \t,r,k) h-> -v^(t,r, -k) 



(9.22) 



uW(t,r,k) i — y u (s) (t,r, -k) 
Combining the two symmetries we see that the action of three-dimensional inversion 
TR : ^(t, x, r) h> ^ (t, x, r) = -T% v i/; v (-t, -x, r)* , (9.23) 
where X = diag{— 1, —1, —1, 1), induces the transformation 



Ti? : 



u (s) (£,r,k) n(- s )(-t,r,k)* 
;W(t,r,k) i/~ s )(-£,r,k)* 



(9.24) 



This is the u^ s ' — > (u^* 1 )* symmetry of the equations (4.15) that we mentioned 
earlier. 

Although TR is a symmetry on the bulk equations we note that this operation 
induces a symmetry of the boundary correlator £y only if the boundary conditions 
also respect the symmetry. The action on the three-dimensional boundary theory can 
be derived from the embedding of the d = 3 Clifford algebra in the d = 4 one via the 
decomposition ( A.4[ ) and the decomposition of the asymptotic expansion variables 



Aa = (0, a a ), Ba = (6 Q ,0) described in section 5.2 Given (9.23) we find 



TR : 



o«(p) ^ -<(P) 
&a (p) ^ &*(p) 



(9.25) 



Hence, given the definition (15.211) and since the basis vector-spinors e*a are real, we 



have the action on the correlation function 



t u (u, k) i-> tn(u, k)" 



(9.26) 



and similarly for £22- 

Let us now examine the action of the TR transformation on the horizon boundary 
condition. We will argue that at zero temperature, as the frequency goes to zero, it 
becomes a symmetry. This limit is subtle because of the irregular singular point at the 
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horizon. However here we can use our small u analysis of the previous section. We saw 
that, for small u, ingoing boundary conditions imply that solution in the AdS 2 x M 2 



region can be written in terms of Whittaker functions as in (9.9). In the overlap 



region, the corresponding near boundary expansion is given in (9.11 ). In this region it 
is then straightforward to take the zero u limit since Qr(u, a(s), v^) ~ (uL^ 2 )) 2uk 
for all s. In the limit we thus have 



x V {r) = (r _ r+ )-*+"> ( M (3/2) Si + U ™ S3 ) + . . . , 



,(-V2). 



(9.27) 



where the constant vector-spinors s\ and S3 are defined in (8.20) and have the form 



rj 



Cie n 



V 



S3 



ry - c 2 e 



rf 



(9.28) 



where c\c* 2 = 1. Note that for the strictly u> = equations of motion, the essen- 
tially singular nature of the horizon is absent and there are four independent solu- 



tions (9.12). What the argument above implies is that at zero frequency we should 



retain only solutions 77(1) and 77(3) that are regular at the horizon (see also JT 57 
(Note that this is not the case for the zero-frequency limit of the T > solutions 
discussed in section [6j) The gauge parameter similarly has two types of solutions 
e ~ (r — r + )~^ ±Uk and, by a similar argument, we discard the lower sign, keeping the 



regular solution. Feeding the expansion through (4.14) to obtain the gauge variations 
of the gravitino we find that we can gauge away either of u^ 2 ^ or Uq l ^ 2 ^ completely. 



Given the action (9.24) we have, under the inversion map, 

1 



STR) 



(TR) 



r] — c x e 

— * H 



7/ 
rf 



c 2 
1 

Cl 



-s 3 



(9.29) 



Hence, in the overlap region, 



x v(TR) {r) 



r — r 1 



u, 



(-1/2). 


Cl 



u, 



(3/2)* 



c 2 



-S3 



+ ... 



(9.30) 



which is the same form as the original solution. Thus the TR symmetry respects 
the ingoing boundary conditions in the zero u limit. This fact relies crucially on the 
conjugation property of ci y2 . More precisely one can show that the TR-transformed 



solution is, up to linear rescalings, gauge equivalent to the original solution (9.27) 



Since the Rarita-Schwinger equation is a linear equation, equivalence up to linear 
rescalings is sufficient to ensure equality of the correlation functions at infinity. The 
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numerator and denominator of (5.26) and (5.27) scale homogeneously with such linear 



factors and so they cancel in the final result. 

We have thus shown that the TR action is a symmetry of the system of equations 
and the horizon boundary conditions at zero temperature, in the limit u = 0. This 
implies the condition 

t n ((U) = t n ((ur (9.3i) 

on the retarded propagator of the boundary theory. An equivalent argument, involv- 
ing the p sector equations can be made for £22(0, k). Thus, at zero frequency, the 
retarded Green's function is purely real and hence the spectral function A(0, k) = 
lm£u(0, k) vanishes. 

In fact we can make a more general statement about the reality of certain expan- 



sion coefficients in the retarded Green's function (9.17). Since the boundary condition 



at higher orders in perturbation theory is fully determined by the boundary condition 



at order zero (consider the expansion of (8.16), or alternatively (9.14), in powers of 



ui£; also see (5j ). we deduce that if the latter is TR invariant, so is the former. Note 
now that the arguments above also show that the quantities m and hf n , at n-th 
order in perturbation theory of the outer region, obeying the boundary conditions 



specified in (9.12) are mapped under TR into 



ctf^(k) 



*?(!>(*) 



(9.32) 



An exactly analogous argument can be made for the latter boundary condition in 



(9.12) involving S2 , S4 implying that 



*t,(2) 



%{2) 



(k) 



(k) 



%(?) 



(k) 



(9.33) 



If we now combine this fact with the result (9.17), we can see the numerically 



observed power- law scaling for low frequencies analytically. Expanding (9.17) for 
small oj we find to leading order 



t u (u, k) = £ n (0, k)ll + J2 C n (k)(co£) n g(u, v k ) + D n (k)(u£) n 



(9.34) 



n=0 



where the C n (k) and D n {k) are fc-dependent constants depending on the (^-independent 



coeficients auybuy To be explicit, expanding the numerator of (9.17) in powers of 



10 



For this computation it is convenenient to make the gauge choice p ■ a = 0. Note that we should 



use the relation (5.24 1 to find the TR action on quantities like a^i). Note also that all quantities 
including the e l a and p 2 have to be expanded in powers of lu£. 
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oo£ gives the first two coefficients 

C x {k) ~- 



™ 6 °,(2) ^1(2) 



D x (k) = — p-_^ + 7 ^, (9.35) 



A) 



where pg = p 2 | w =o and the action of TR (9.25) induces the simple map 

Ci(Jfe) C x {kY , D t {k) D^ife)* . (9.36) 

Since we saw above that at = the boundary conditions are TR invariant, we can 
conclude that these coefficients are real. In fact, a little more thought shows that 
all higher coefficients D n (k) and C n (k) are all mapped to their complex conjugates 
under the action of TR. Hence we conclude that they also are real. Finally, putting 
these arguments together with £n(0, k) G R, we can extract the scaling relation at 
leading order 

Im(* u (w,Jfe)) oc tu 2uk , (9.37) 

valid for the spectral function near the origin, exactly as we saw in our numerical 
results. 



10 Quasi-normal modes and their physics 

In this section we establish the location in the u, k plane of the quasi-normal modes 
(QNMs) of the gravitino in the AdS-RN black-brane geometry. Recall (e.g. |7i~] ) 
that the QNMs should have ingoing boundary conditions at the future horizon and 
normalisable boundary conditions at the AdS boundary. The latter condition implies 
that the QNMs are associated with poles of the Green's function when considered 
as functions of complex u. Furthermore, the QNMs with u close to the real axis 
are associated with prominent features of the Green's function when considered as 
functions of real u. Determining the QNM spectrum is thus a non-trivial check on 
the structure of the supercurrent two-point function we studied above. The results 
of this section, to our knowledge, are the first determination of the QNM spectrum 
of the gravitino in an AdS black-brane background. Our approach is numerical and 



makes use of Leaver's method as first developed in 72 and more recently applied to 



AdS/CMTp^jin (9 , 73 . For simplicity we will just deal with T ^ 0. Our results here 
are an independent confirmation of the pole structure in Figures fTMJ 



11 We thank the authors of j9 73 for sharing technical details of their work with us. 
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10.1 Quasinormal modes for T > 



a) 



c) 



0.0- 
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Figure 5: i/ere we s/iow t/ie highest -lying quasinormal frequencies in the complex u 
plane. Background parameters are chosen so that temperature T/fi = 0.08. We set 
Ik = 2.1 in panel a), Ik = 1.7 in panel b), £k = 1.1 in panel c) and Ik = 0.1 in panel 
d). One can clearly see the pole moving to oj = as k — >■ 0, corresponding to the 
supersymmetric sound pole at finite T. There is a line of poles just off the imaginary 
axis, which becomes a branch cut in the zero-temperature limit. We note that the 
Leaver method does not give any information about the strength of the pole. The 
figures show density plots of |detM'(ix>)/detM((X>)| ; where the matrix M(u) is defined 
in MJfy. 



We now turn to describing the appropriate ingoing boundary conditions at the 
horizon. At finite temperature the horizon at r = r + is a regular singular point of the 
gravitino equations and we can use the expansion studied in section [6] above. In our 
implementation of the Leaver method we shall specifically not choose a gauge at the 
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horizon (having fixed the gauge near the boundary), which is the opposite approach 
to that employed in the calculation of the spectral density above, where we fixed 
the gauge at the horizon, but not near the boundary. These are merely matters of 
convenience, of course since the physical quantities we extract, such as the location 
of the QNMs, are gauge-invariant objects. 




a) Rete#3i b) ReHu) 

Figure 6: Here we show the highest-lying quasinormal frequencies in the complex u 
plane for kl < 0. Background parameters are chosen so that temperature Tj \x = 0.08. 
We set Ik = —1.5 in panel a) and Ik = —2.8 in panel b). We see that the phonino 
pole bounces off the real axis at kl = and moves back down into the negative Im u 
plane with Re a; < — \x. This corresponds to the spread of the high spectral density 
region associated with supersound to negative u at negative ki in Fig. [#| 

We now have assembled all the information necessary to set up a version of 
Leaver's method, adapted to our system. In this approach one factors out the de- 
sired leading behaviour at the horizon as well as infinity and expands the solutions 
around the midpoint of the interval between horizon and infinity. For this purpose it 
is very useful to change coordinates to Z = £ 2 /r and we shall do so in this section. In 
addition, we find it convenient for the numerics to set £ — r+ = 1, using the scaling 
symmetries of the equations, so that the temperature is entirely set by the chemical 
potential \i. Thus, in detail, we expand our vector spinor (in the rj sector) to some 
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arbitrary, but fixed order N as 

N 

u (3/2) = ( Z -l)-I-fz 3 / 2 ^ai 3 / 2 )(Z-|r, 

n=0 

TV 

= {Z-l)-\-^Z*l 2 Y, a n /2 \ Z -\Y, 

n=0 
N 

= {Z-l)^Z*l 2 Y, a{ n m ( Z -\Ti 
n=0 
N 

u (-3/2) = ( Z -l)f-fz 3 / 2 ^a(- 3 / 2 )(Z-ir, . (10.1) 

n=0 

We then substitute this ansatz into the gravitino equations of motion and collect the 
resulting linear equations into a matrix equation of the form 

H .M 3 / 2 1 

£ £ M -(»,q,T)^') = 0, (io.2) 

s '=-3/2m=0 

for complex frequencies tt> where, as before, to = u/27rT and q = k/2nT '. Solutions 
of this equation are forced to be trivial, unless the condition 

= 0, (10.3) 

is met for some special complex frequency tt>*. The determinant is taken over all 
indices including the helicity index s. Solutions to this equation are the quasinormal 
frequencies of the background, under spin-3/2 perturbations. In practise we extract 
the matrix using computer algebra and proceed to numerically solve the deter- 
minant condition. 

We present our results for a finite temperature T/fi ~ 0.08 in Fig. [5j The 
behaviour of the highest QNMs at other finite temperatures is qualitatively very 
similar and we do not present explicit plots here. We see that there are two highest- 
lying poles, one in the negative u plane and one in the positive u plane. The former 
corresponds to the phonino pole and approaches the real axis as k — > and Co — > 0, 
while the latter corresponds to the second region of high spectral weight identified in 
Fig. [2] We see clearly how, as k is increased, these two poles cross over so that for 
lower momenta the phonino pole dominates, while for low k it is the latter pole that 
resides closest to the real axis, consistent with it dominating the spectrum. We also 
see another set of poles going down the complex u plane near the axis (but just off 
it). It should be noted that the representation of the poles as computed in the Leaver 



det 



Km (»*,q,T) 
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method can be slightly misleading as the strength of the poles on the plot is not 
related to its residue, as it would be, for instance, if one computed the full spectral 
function in the complex u> plane. The poles going down the imaginary axis are in fact 
much weaker than the two leading ones we discuss above and correspondingly do not 
give an appreciable imprint on the physical spectral function of the supercurrent. A 



similar point was made, for the case of a bosonic field in 74 . We also studied the 
QNM structure at T = 0, and found results that are consistent with the corresponding 
spectral density presented in Fig. |4| 

11 Final comments 

In this paper and its companion [I] we have carried out a detailed analysis of the 
spectral function of the supercurrent for the general class of d = 3 , N = 2 SCFTs 
which have D = 10 or D = 11 supergravity duals. This was achieved by using the 
consistent KK truncation results of [58] and then solving the gravitino equations of 
iV = 2 D = 4 minimal gauged supergravity in the AdS-RN black-brane background. 
A major feature of the spectral function is the presence of the phonino pole at k = 
and u + n = 0. We calculated the dispersion relation for this pole deducing a speed 



consistent with the hydrodynamical results of 65 and also the diffusion constant. It 



would be interesting to go beyond our numerical results for the dispersion relation and 



obtain its value analytically. It would also be interesting to extend the analysis of 65 
to finite chemical potential. A second major feature of the spectral function is the 
depletion of spectral weight around u = 0, which is accentuated for low temperatures. 
Indeed at zero temperature we showed that the spectral function vanishes at oj = 
giving rise to a power-law soft gap that is controlled by the locally quantum critical 
point dual to the AdS2 xR 2 solution. In the latter background we obtained an analytic 
expression for the Green's function by solving the gravitino equations exactly in terms 
of Whittaker functions. 

It is noteworthy that there is no Fermi surface in the supercurrent correlator 



at (bulk) tree level, a result that has been recently confirmed in 57 . While it is 
possible that this indicates that there is no Fermi surface in these top-down models, 
we think it is more likely to be a consequence of the particular correlator being 
considered and in particular a Fermi surface will be seen in other Fermion correlators 
and/or at higher loop level. We hope to analyse a concrete top down example in the 
future to develop a better understanding of this issue. One possibility is to study 



the top-down models given in 75 and furthermore examine the behaviour of the 
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supercurrent and other fermion correlators in the superfluid phase that appears at 



low temperatures 76 77 . We saw in section 8 that the conformal dimensions of the 
operators in the one- dimensional conformal field theory in the far IR, dual to the 
AdS 2 region, are all real. This means that the "log-oscillatory" behaviour seen in the 
bottom-up models of [3j[5] is absent. In the bottom-up models this behaviour was 
correlated with the presence of a Fermi surface and one might wonder if there is some 
deeper connection between these two phenomena. It would be worthwhile extending 
our calculations to the alternative quantisation of the gravitino field to check which 
of these features persist. Although this just entails changing the AdS^ boundary 
conditions, some care will be required since the alternative quantisation breaks the 
boundary supersymmetry and hence the gravitino will have additional propagating 
modes in the bulk. 
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A Spinor conventions 

Throughout we have chosen a mostly minus signature convention. For the d = 3 
gamma-matrices 7° we will use the explicit real basis 

7* = ia 2 , 7 X = cxi, 7^ = a 3 , (A.l) 

with , y txy = +1. For a d = 3 spinor s we define s = s'7*. Parity in d = 3 will be 
taken to act on the coordinates via (t, x,y) — > (t,x, —y) and acting on spinors via 

s(t,x,y) -+7 y s(t,x,-y). (A.2) 

With this definition observe that the bi-linear s^s transforms as a vector. However, 
the bi-linear ss transforms as a pseudo-scalar (a scalar bi-linear can only be con- 



structed from a spinor transforming as in (A.2) with another spinor that is chosen 
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to transform under parity with an extra overall minus sign). We will also define the 
basis of d = 3 spinors 

m =(j), »=(;), (a.3) 

which satisfy / ~f tx m = m, / -f tx n = —n and fhn = —nm = 1. 

For the .0 = 4 gamma-matrices, with tangent space indices, we will choose the 
real representation 

r d = 7 Q ®(7i, r = l<g><7 3 . (A.4) 

In the four-dimensional space spanned by coordinates (t, r, x, y) we have projectors 
made from the following chirality operators 

p(2) = T i T r p(3) = 1 p(2)^ = T if & , A5 s 

|k| 

Since these commute, [r^ 2 ) , T^ 3 )] = 0, we may label the components of four-dimensional 
spinors (and vector spinors) in terms of simultaneous eigenvalues under them. Let us 
denote their simultaneous eigenspinors as 

T(V = ^ r(V = "P ± , T^r ] ± = ±r ] ± , T^p ± = ±p ± . (A.6) 
We find it convenient to work in a basis with 

\ m I \ m I \ n I \ n 

We also have 

rV 

r V : 
rV : 

rV 

B Positivity of the Spectral function 

By definition the retarded correlation function in position space has the form 

G a p(x -y) = i6{x° - y°) ({S a (x), S (y)}) . (B.l) 

If we define 





rV = ± p t , 




-if , 


r V = -p t , 




± v ± , 


r V = TP ± , 




TP±, 


r V = tv* ■ 


(A.8) 
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then, after a Fourier transform, the spectral function is given by 

A a/3 (uj,k) = lm.G al3 (u, k) = r]'^(u,k), (B.3) 

where r}'^(oj, k) is the Fourier transform of fj'^(x — y). 

Formally we can write a complete set of states of the Hilbert space as 

1 = J dti(q)dti(X)\q,X)(q,X\ , (B.4) 

where q = (_E q , q), with = ^/q 2 + m(A) 2 , labels the total three-momentum and A 
is a formal label for the other degrees of freedom, while dfi(q) = d 2 q/2E (l (27r) 2 and 
dfi(X) are the appropriate measures. Taking the thermal expectation value at finite 
chemical potential /i, we have, by inserting a complete set of states and tracing, 

(S a (x)S p (y)) = J dpi{q)dpi{X)d^q')d^X')e-^ + ^ N 

(q, X\ S a (x) \q', A') (q', X'\ Sp(y) \q, A) 

(B.5) 

d^(q)dfi{X)dfi{ q yfi{Xy- pE+ ^ N e- iiq - q ' )ix - y) 
(q,X\S a (0) \q',X'} (4,\'\Sp(p)\q,\), 
where E = q° is the energy of the state \q, A) and N is the .R-charge. We define the 



matrix element that appears in (B.5) 



M a (q, q\ A, A') = (q, X\ S a (0) \q\ A') . (B.6) 

The gamma-tracelessness and weakly gauged conservation of the supercurrent S a 
implies that 

j a M a = 0, p a M a = 0, (B.7) 
where p = (a), k) = (q° — q'° + /i, q — q'). Thus we can expand M a in terms of the 



basis (2.7) 



M a (q, q>, A, A') = Zi{q, q', A, A')e«. (B.8) 

Changing variables to the three-momentum q — q' so that the measure can be written 

as 

diM(q)diMtf)dii(\)dfi(\') = dfi(q, q', A, A') - , (B.9) 

we have 

€^ -V) = J d -^0 1 Mp)efef e ~^ ){X ~ y) > (B.10) 



49 



where 

MP)= I dj,(q,q',\,\')e-^ N 



(B.H) 

x [Z i (q,q',X,X')Z*(q,q',X,X') + Z i (q',q,X',X)Z*(q',q,X',X)] . 



In the notation of section [2] we have Aij(p) = Imtjj(p). In particular we deduce that 

A(p) =lmtu (p) =A u (p)>0, (B. 12) 



because the integrand in (B.ll) is positive definite. 



C The Rarita-Schwinger equation in AdSd+i 

In this appendix we consider the Rarita-Schwinger equation for an uncharged grav- 
itino in d + 1 spacetime dimensions. We found this to be a useful exercise in order to 
elucidate some of the analogous calculations for the charged gravitino in 3 + 1 space- 
time dimensions. The material presented here may also have other applications. 

The Rarita-Schwinger equation for an uncharged gravitino in d + 1 spacetime 
dimensions is given by 

777 

r^v^ p --r^ = o, (c.i) 

where V is the Levi-Civita connection. At this stage the parameter I > is arbitrary, 
but is inserted for convenience, since shortly we will restrict to the case of AdSd+i 
with radius taken to be £. For definiteness we will take m > 0. The case that is most 
relevant to the bulk of this paper is d = 3 and rh = —tm = +1. 



Following 78 we write 
whence upon contraction, first with V M , then with T^, we obtain 

771 

o = yy(r.7/;)-v 2 (r-7/;)-y(v-7/;) + v^yr)--(yr-7/;-v-7A) 

777 n 

o = (d-i)[y(r^)-v-v]-^pr-v>. (c.3) 

We can use various curvature identities to remove the double derivative terms from 
the first expression. Since t/> m is in the spin-3/2 representation we have 

[v„ v u ] r = R^ pa ^ + \R^ 5 T aS r , (cu) 
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from which we can deduce 



yy(r-vo = v 2 (r-^) - \r{v-^) 



(C.5) 



As advertised, these identities allow us to eliminate all two-derivative terms from the 



first equation in (C.3) 



The resulting expressions are particularly simple in AdSd+i spacetime, for which 
R^ v = —(d/ l 2 )g^ u . Indeed after a bit of algebra we conclude 



0= ( ; ,r- (d V ) (r-V). 



(C.6) 



Thus, we have to distinguish two cases. If the gravitino is physically massless, i.e. if 

(d-1) 



m 



this equation is satisfied identically. If not, we must set 



r-^ = o. 



v- v = o 



(C.7) 



(C.8) 



where the second constraint follows from the second equation in (C.3) once T -ip = 



is imposed. If the gravitino is physically massless, as in the supersymmetric case we 
consider in the bulk of the paper, we can in fact use the local supersymmetry, to 
impose these same two constraints, taking note, however, that they do not fix the 
gauge freedom completely. More specifically, the supersymmetry transformations are 
(setting all other bosonic fields in the supergravity to zero) 



(C.9) 



c%= (£> A + -r A le, 
and the residual supersymmetry transformations read 



m + 1 , 

# + —7— |e = 0, 



(CIO) 



The equations Eq. (C.8) project out the spin-1/2 parts of the Rarita-Schwinger 



equation and we thus see that solving (C.l) is equivalent to solving 



v-v = o. 



(C.11) 
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C.l Asymptotic analysis of Rarita-Schwinger in AdSd+i 

Let us write down the asymptotic form of these equations at the boundary of AdSd+i- 
This allows us to deduce the leading asymptotic behaviour. We write the AdSd+i 
metric in Poincare coordinates 

ds 2 = ^+^dr 2 , (C.12) 

where ds d is the flat metric of The equations can be succinctly written as 

r (rd r + \d) ipjx + jt)^ - m% = -T^f , (C.13) 

where D a is the connection on ds 2 ,. If we have a power series of the form ip^ = 
^2 p C^r~ p the equations imply 

[(p - Id) r - m] = rtdt x) + r*c« 

[( p _I rf _l)r-m]^ p) =£ 2 ^ p - 1) . 

We now have to distinguish two cases, the first of which is more relevant for this 
paper. 

fa £Z+ | 

Then we can develop an asymptotic series the leading behaviour of which is 

^ = r A " d 4 0) (1 + 0(r-i)) + r^Bf (l + 0(r^)) , 
^ = r A ~ d - l AP (1 + ^(r- 1 )) + T-^BP (1 + ^(r- 1 )) , 

where A = \d + fa and 

(u+n4 0) = o, (i-n4 0) = o, 

- m - m (C.16J 

(u + r r )4 ) = » (u-r r )s« = o. 



Note from (C.14) that successive terms have opposite T r chirality. Thus the B®) 



series does not mix with the A^' series unless fa e Z + |. 

In this case B<-°) and A^ 2 ™) have the same chirality. To get the full solution we need 
to include a logarithmic term so the asymptotic solution takes the form 

^ = r -A + 2^ A (0) ^ + Q{r -X^ + r -A lnri (0) ^ + ^ 

^ = r -A+2 A -l A (D ^ + 0(r -l^ + r -A-l lnr ^(D p + 0(r -l^ ^ ( C ' ' ' ' 
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(rd r + d+±)fa + -D-iP = 0. (C.20) 



where 

(u+n4 0) = o, (i-nif = o, 

(i + r)4 =o, (iL-r r )4 } = o. 

Next we consider the conditions F ■ ip — and D -ip — 0, focusing on the m ^ Z+ 1 
case. To leading order the gamma-traceless condition is simply 

r & A (o) = T&B (o) = Q ^ ^ c ig ^ 

that is gamma-tracelessness on the boundary. The D -ip — condition reads, after 
using the r ■ ip = condition, 

r 

To leading order this implies 

(l(d-l)+m) Af ] = -t 2 D ■ m , 
(\{d-l)-m)Bf ] = -l 2 D-B^ . 

Thus provided m ^ \{d — 1) we see that (gamma-traceless) 4 and i?^ are the 
independent degrees of freedom. 

The gravitino is dual to an operator with dimension A in the dual CFT. The 
4°^ and coefficients specify the deformations of the CFT by sources and the 
expectation values of the dual operator, respectively. For the rh G Z + \ case the 
coefficient specifying the expectation value behaviour is just the 2m-th term of the 
(non-log) series in (C.17) i.e. we can identify B^ = A^ m \ 



Finally we consider the supersymmetric case where rh — \{d — 1), again focusing 
on the rh ^ Z + | case. First note that the asymptotic expansion now reads 



fa = r^Af (1 + Oir- 1 )) + r- d+1 / 2 4 0) (l + <D(r~ 1 )) 
fa = r- 3 / 2 4 1} (1 + Oir' 1 )) + r-^B^ (1 + Oir- 1 )) 



(C.22) 



The gamma-traceless condition again implies 

r&A (o) = T a B (o) = _ ^ C 23 -j 

The D ■ ip = condition implies 

D-B {0) = 0, (C.24) 
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and is now undetermined. To elucidate what is going on we need to consider the 
residual supersymmetry transformations. The residual supersymmetry parameter e 
satisfies (CIO) which becomes 

r (rd r + \d)e + \{d + l)e + —be = . (C.25) 
The asymptotic solution is then 

e = r VV-D (1 + Oir- 1 )) + (i + ( r -i)) ? ( C . 2 6) 

where 

(ii+rV 1) = 0) (u-r f )6 (1) = o. (c.27) 

We also have 

a (o) = _l_Z) a (-i). (C .28) 



The supersymmetry variation of ^ in (C.9) is given by 

S^ & =^T & (l + r)e+^D & e 

r 1 

St/ft = jd r e + — T f e . 



(C.29) 



We find the leading variations are (note that the potential SA f , SA& ^ and SAf^ 
terms all vanish) 

8A® =t(f)*- -T A £>) a™ 



d a 



(C.30) 



SB® = -rV6« . 

We find that A? is gauge-dependent (though remains gamma-traceless) as is -B^. 
We can always then choose a gauge where D ■ A^ = and B^ = 0. The final result 
matches with our expectations for the supersymmetric case. The coefficient giving 
the expectation value, BjP \ is gauge-invariant. On the other hand the coefficient 
fixing the source, A% , is gauge- dependent, corresponding to a weak gauging of the 
supersymmetry in the dual SCFT. 



D Comments on our numerical approach 

Here we describe briefly some aspects of our numerics, which are designed so as to cope 
with the highly oscillatory behaviour of the equations, especially at low temperatures 
(including T = 0). 



54 



The general idea is to integrate equations (4.15) from the horizon at r = r + out 
to infinity (or, in practise, to a very large value of the radial variable) and then to 
extract the boundary quantities and b& of Eq. (5.20) by fitting the resulting data 



to the asymptotic expansion of section [5j However, especially for small T / //, the 
near-horizon behaviour of the gravitino components is highly oscillatory, aggravated 



by a different singularity structure (see Eqs. (6.7) and (7.4)) at the horizon for 
different helicity components and thus different rates of oscillations and divergence. 
This means that the equations are often stiff in the numerical sense and as a result 
once one has integrated a long way out towards infinity the numerical error has 
accumulated to an extent that the analytical expansions of section [5] are no longer a 
good fit. Since we need to extract the leading order coefficient at 0(r~ 3 / 2 ) as well as 
the expectation value at 0(r~ 7 ^ 2 ), down by two orders in the expansion, as well as 
the intermediate orders (in order to verify gauge invariance) a more reliable method 
is needed that works to very high accuracy. 

We circumvent the problems outlined above by side-stepping direct integration 
in favour of a shooting algorithm. We now briefly summarise our procedure, which 
works at T = and T > 0. We present the procedure for both r\ and p sectors at the 
same time, but in practise, since they completely decouple in the equations (4.15) it 
is easier to just focus on one at a time. Whenever we do numerical calculations we 
use scaling symmetries to set I = 1 and r + = 1 so that the effective temperature of 
the background is fixed by the sole remaining dimensionful parameter /i. We also 
assume that to ^ for the steps below. 

1. Use computer algebra to develop the near-horizon expansion to a very high 
order. Both at T = and T > we choose ingoing boundary conditions 

jf, ~ e -^*-^ ; (D.l) 

using the appropriate tortoise variable r* for each case. We explicitly implement 
the gauge transformations on the series so that we can generate from a given 
choice all other solutions which are in the same gauge orbit. One check on our 
numerics is, then, that all gauge-equivalent near-horizon boundary conditions 
give the same physical correlation functions at infinity. There are four complex 
numbers specifying the general near-horizon solution, which we can interpret 
as specifying initial conditions for the two physical helicity states as well as the 
gauge parameters in the rj and p sectors. 

2. Use computer algebra to develop the near-boundary expansion of the gravitino 
to a very high order. This expansion is specified by eight complex numbers, 
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comprising the independent degrees of freedom in the vector spinors and 
four in the rj sector and four in the p sector. 

3. Integrate in from infinity and out from the horizon and match the eight inde- 
pendent components of the gravitino vS s \v^ , s = — §, — |, |, | at the midpoint, 
which is of course arbitrary, but can be chosen at a numerically convenient lo- 
cation. We find it numerically more favourable to use a variable Z = —£ 2 /r 
to do the actual integrations, whence the horizon is at Z = 1 and asymptotic 
infinity is at Z = 0, so that a convenient and obvious choice for the midpoint 
is Z mid = 0.5. 

4. We then use a shooting algorithm to match the solution, z m id, which is equiva- 
lent to determining the data at infinity as a function of the near-horizon data 
as well as the single parameter \x determining the effective temperature of the 
background. At every iteration step of the shooting method we must numeri- 
cally invert an 8 x 8 matrix, which can be easily and stably done numerically. 

5. Run the above steps for a fine grid of values of u (if desired complex) and k 
and store, e.g. as binary files. We can subsequently use the stored collective 
data to assemble the various plots and fits shown in the bulk of this paper and 
the companion paper [1]. The data for the correlators presented in these works 
is based on (9(500, 000) data points, each obtained by the shooting algorithm 
above, which converges after a few steps to an accuracy of O(10~ 10 ). This gave 
a total running time on a top-range commercial dual-core processor machine 
of approximately 1 week (in practise several machines with several cores were 
used in parallel). 
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E Asymptotics of Whittaker Functions 



For small values of uz we have 



-2iuzf^ F( , 2 " fc) , (1 + CM) 



+ (-Ma*)*"* wl ~ty~ \ (1 + OW) 



2za;L (2) ) 1 / 2 -^ 



T(2z/ fc ) 



r (i - a(s) + i/ fc ) 



L (2) 



;i + 0M) 



+ h— ^(a(s),i/ fc )(l + 0(a;z)) 

L (2) 



(E.1) 



where, in writing the expression in the equality in the fifth row, we could hide the 
leading Euler gamma functions by choosing a convenient overall normalisation, as we 
do in the bulk of the paper. The ratio of gamma functions Q R (a(s), u k ) is the retarded 
correlation function of the one-dimensional Spin(l, 1) representation labeled by s: 



T(-2u k )T (l - a(s) + u k ) 

g R (a(s), u k ) = A- k) 2 \> k M2uL [2) r* . (E.2) 

T(2u k )T (2 - a(s) - u k ) 



In a similar manner, by considering the asymptotics of the outgoing Whittaker func- 
tions W- a ( s \ Uk (2iuz) we can derive the advanced correlation function for each helicity 
label s: 

GaHs), u k ) = e™ fc ~ \ ^' 77 ' ~: j (2uL m y»* . (E.3) 



T{2u k )Y{\ + a{s)-u k ) 



It is interesting to note that 

Q R (u, a(s), u k ) e~ 2 ^ - e -^ L (2) 



(E.4) 



Q A (UJ a(s),U k ) e 2i-KU k _ e -V2nqL (2) 

where we have used the fact that for all of the helicity states of the gravitino we 
have e 2ms = — 1, a reflection that these are fermionic states. Quotients of Q for two 
different helicities s 7^ s' are similar and can be related exactly to the above result 
by using the recursion identity for the Euler gamma function. 
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